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ABSTEACT 



The external nucleate pool-boiling heat- transfer coeffi- 
cient of a horizontal smooth copper tube in R- 11 4-oil 
mixtures (0 to 10 percent oil) was measured for heat fluxes 
from 1 to 100 kw/m^ at two different saturation temperatures 
(-2.2 °C and 6.7 OC) . A copper-nickel tube coated with the 
Onion Carbide "High Flux" coating was similarly tested. The 
High Flux coating was found to improve the heat-transfer 
coefficient by at least a factor of 7 in oil- free R- 1 1 4 . Oil 
resulted in about a 20 percent reduction of the heat- 
transfer coefficient of the High Flux surface at heat fluxes 
less than 30 kW/m^ and up to an 80 percent reduction at heat 
fluxes above 30 kW/m^ with greater than 6 percent oil. 
Under all conditions, the High Flux coated tube outperformed 
the smooth copper tube. 
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I. IHTBODOCTION 



a. BACKGEOOHD 

The U.S. Navy currently uses refrigerant R- 1 1 4 in 
centrifugal chilled-water air-conditioning plants aboard 
submarines and surface ships. The Navy hopes to reduce the 
size of these units and increase their performance by using 
enhanced evaporator and condenser surfaces. An experiment 
by Arai et al. [Ref. 1] produced a prototype 200- Ton R-12 
centrifugal water, chiller that was 28 percent shorter in 
length and had a 50-70 percent improvement in the overall 
heat- transfer coefficient by employing the enhanced surface 
’’Thermoexcel E” made by the Hitachi Company. Comparisons of 
various enhanced commercial tubes by filmaz and Westwater 
[Ref. 2'], Marto and Lepere [Ref- 3], and Carnavos [Ref. 4] 
for various refrigerants other than R-114 indicated that a 
porous-coating-enhanced surface, such as Union Carbide’s 
"High Flux,” will exhibit the best boiling heat- transfer 
performance in a pure refrigerant. 

The High Flux surface (see Figure 1.1) consists of a 
sintered metallic matrix bonded to a metallic substrate. 
The surface is produced by coating a smooth tube with a 
binder-solvent mixture and then applying a mixture of metal 
and braze alloy powder; the tube is placed in a furnace to 
evaporate the solvent, binder, and melt braze alloy thus 

forming a porous structure having multiple reentrant cavi- 
ties to enhance nucleation. [ Ref. 5 ] 

Since the High Flux surface will be employed in a 
refrigeration unit using an oil-lubricated, hermetically- 
sealed, compressor, some amount of oil is always present in 
the evaporator. Studies by Henrici and Hesse [Bef. 6] for 
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Fi gu re 1.1. 



Schematic and Scanning Electron Micrograph 
(500x) of High Flux Surface. 



smooth tubes and Stephan £Ref. 7] for the Gewa-T surface 
(manufactured by the Hieland Company) in E-1 14-oil mixtures 
indicate that the heat-transf er coefficient of enhanced 
surfaces can be significantly altered by oil. 

Experimental data showing the effect of oil concentra- 
tion on the heat-transfer coefficient of the High Flux 
surface in E-1 14 are lacking, thus motivating the present 
investigation. This investigation was funded by the David 
W. Taylor Naval Ship Research and Development Center. 
Details of the experimental apparatus are described by 
Karasabun [Ref. 8]. The smooth copper tubes were supplied 
by the Hieland Company. The High Flux coated copper-nickel 
tubes were supplied by the Onion Carbide Corporation. 

B. THESIS OBJECTIVES 

The objectives of this thesis are: 

1. Take boiling data on a smooth tube in R-114 with and 
without oil for comparison with the data of other 
researchers, and to provide baseline data for evalu- 
ating the boiling performance of the High Flux tube. 

2. Take boiling data on a High Flux tube for various oil 
concentra tions (0 to 10 percent by mass) . 

3. Study the effect of saturation temperature on the 
R-114 boiling behavior. 

4. Attempt to sample oil locally in the near vicinity of 
a tube to investigate the possibility of an oil 
concentra tion gradient around the tube during 
operation. 
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II. EEVIEH OF BE FHIGE RA NT-OIL HIXTD^ BEHAVlOB 

A. HOCLEATE BOILIHG OF BEFfilGEEANT-OIL MIXTOBES FROM SMOOTH 

TUBES 

In 1963, Stephan [fief- 9] published a milestone paper on 
the influence of oil on the boiling heat transfer of R-12. 
The effects he noted have been observed in most refrige- 
rants, including R-114. In 1972, Henrici and Hesse [Ref. 6] 
updated Stephan's work for R-1 14-oil mixtures boiling from a 
smooth copper tube. Figures 2-1 and 2.2 summarize Henrici 
and Hesse's results. Figure 2.1 shows that oil generally 
lowers the heat- transfer coefficient, and that at high heat 
fluxes and high oil concentrations (10 percent), the effect 
grows more pronounced (slope decreases) . Figure 2.2 shows 
that at some oil/heat flux combinations, the heat-transfer 
coefficient may actually be improved by the addition of oil. 
Chongrungreong and Sauer [fief. 10] suggest that the heat- 
transfer behavior of refrigerant-oil mixtures can be attrib- 
uted to 5 major factors; 1) the physical properties of the 
refrigerant-oil mixture, 2) the saturation temperature (or 
boiling pressure) , 3) the tube diameter, 4) the surface 

condition of the tube (roughness) , and 5) the hydrostatic 
liquid head above the tube. 

T • Physica l Prope rties 

Refrigerant- oil mixtures have significantly 
different physical properties than pure refrigerants. 
Jensen and Jackman [fief. 11] report that density and 

specific heat behave ideally in refrigerant-oil mixtures, 
but that viscosity and surface tension do not. Ideal 
behavior of refrigerant-oil mixture density and specific 
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Figure 2.1. Effect of Oil on Boiling Coefficient 
of R-114 (from Ref. 6). 
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heat does not mean linear behavior. The governing egua- 
tions are; 
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used in their 1971 experiment. As shown in Figure 2.3, the 
surface tension of the mixture first decreased up to an oil 
concentration of 2.5 percent, and then increased continu- 
ously with increasing oil concentration. This type of non- 
linear behavior makes explaining the change in heat- tra nsfer 
coefficient of ref rigeraut- oil mixtures, due to the changing 
physical properties of these mixtures, both difficult, and 
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possibly non-general. The behavior of refrigerant-cil 
mixtures is specific to the particular mixture components, 
and may be dependent upon the kind of oil being used. Seme 
qualitative consequences of adding oil to a refrigerant, 
however, can be noted. 

The most observable result of adding oil to refrige- 
rants is foaming. Oil concentrations above 1 percent result 
in significant amounts of foaming from nucleate boiling. 
The foam bubbles form because the R-114 in the R-114-oil 
mixture is more volatile than oil and vaporizes first, 
creating a gas bubble surrounded by an oil-rich layer (see 
Figure 2.4). Since the bubbles are coated with oil film 
with a higher surface tension than the bulk liquid and a 
have lower density, they rise to the top of the liquid. 
Because of their surface tension, the bubbles build up on 
the liquid surface to produce a foam layer. 

This foaming action, which is most pronounced 
between 1-10 percent oil concentration [Ref. 6], may affect 

the heat transfer of tube bundles significantly. For single 
tubes, it is the oil concentration gradient which would seem 
to play the major role, since the foam rises away from the 
tube surface and could only interact with the tube as it 
sweeps by it from the bottom to the top of the tube. 

The general decrease in the heat-transfer coeffi- 
cient upon adding oil to pure refrigerants (recall Figure 
2.1) is subject to many explanations. Thome [Ref, 14], in 
an extensive review cf the literature, reports that the 
first explanation for the decrease in the heat- transfer 
coefficient of mixtures was presented by Van 5?ijk et al . in 
1956. The effect was explained as being the result of the 
evaporation of the more- vola tile components, leaving an 
oil-rich layer with a higher local boiling point, which 
increases the amount of superheat required to continue 
vaporization and bubble growth, thus reducing the 
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heat-transf er coefficient. Stephan and Preusser [Bef. 12] 
demonstrated conclusively that the work of formation of 
bubbles in a mixture is greater than in an equivalent pure 
fluid. They concluded that the mixture heat-transfer coef- 
ficient is lower than for the equivalent pure fluid because 
of the resulting decrease in the bubble population. Any 
variation from a decreasing mixture heat-transfer coefficent 
(like Henrici and Hesse show in Figure 2.2) is attributed by 
Stephan [Bef. 13] to the non-linear variation of the phys- 
ical properties. Stephan proposes that the plot of surface 
tension (Figure 2.3), along with thermal-property varia- 
tions, accounts for the anomalous rise in the heat-transf er 
coefficient between 3-6 percent oil concentration. 
Chongrongeong and Sauer [Bef. 10] state that it is the rate 
of heat diffusion, governed by the thermal properties of the 
oil-rich layer, that limits the bubble growth and that the 
surf ace-tension effects are neglible. 

Thome [Bef. 14] proposes that all of the above 
factors, as well as the viscosity variation, are important 
in explaining the rise in the heat- transfer coefficient for 
some refrigerant-oil mixture and heat-flux combinations. 

In summary, all researchers agree that the physical 
and thermal properties of a refrigerant-oil mixture are 
important factors in explaining the heat-transfer behavior 
of mixtures. 

2 • Saturation Te mperatu re 

It has long been noted that increased saturation 
temperature (i.e., increased boiling pressure) increases the 
boiling heat-transfer coefficient of surfaces in 
refrigerant-oil mixtures. In 1963, Stephan [ Bef « 9] found 
that at high oil concentrations, the heat-transfer coeffi- 
cient of refrigerant-oil mixtures becomes constant with 
respect to the saturation temperature. Stephan proposed 
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that this is because the addition of oil to a refrigerant 
introduces a large diffusion resistance, and that since the 
velocity of diffusion is almost independent of temperature, 
so should the heat- transfer coefficient become independent 
of temperature at high oil concentration. Figure 2.5 shows 
Henrici and Hesse's data on tlie effect of oil and saturation 
temperature in R-114. With no oil, the effect of raising 
the saturation temperature is seen to be a rise in the heat- 
transfer coefficient. With oil, raising the saturation 
temperature is seen to cause a slight drop in the heat- 
transfer coefficient. This effect has not been explained 
yet . 



3 • Dia me ter 



Cornwell, Schuller, and Einarsson [fief, 
that for smooth tube diameters from 6 mm to 30 
nucleate pool boiling heat-transfer coefficient 
refrigerants falls with increasing diameter. The 
tube diameter was correlated by: 
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Figure 2« 5 Effect of Saturation Temperature 
on fi-1 14-Oil Mixtures {from Bef. 6). 



Since equation. (2.-3) depends on physical parameters, the 
effect of oil could he significant. No known measurements 
of the effect of tuhe diameter in refrigerant-oil mixtures 
have, as yet, been performed. 

Surfac e Cond i tion 

Several researchers [Hef. 16 , 17 , 18 ] have investi- 

gated the effects of surface roughness on the heat transfer 
of pure refrigerants. As surface roughness increases, the 
heat-transf er coefficient was found to increase due to 
increased nucleation. Nishikawa [Ref. 18 ] reported the 
effect of a variety of surface roughnesses in pure E-114 
over a range of pressures from 0.294 NPa (42 psi) to 
2.94 MPa (420 psi) . 

5 • Hydros t a tic Effect 

The liquid cclumn above the boiling surface may 
generate large static pressures which will increase the 
boiling point. For E-114 at 0 <>C (32 op), a 0.3 m (1 ft) 
liquid head will raise the saturation temperature about 5 °C 
(9 OF). In large machines, this may be a significant 
effect. However, for small experimental apparatuses, the 
effect is negligible. 

B. HDCIEATE BOILING OF HEFEIGERANT-OIL MIXTDEES FROM 

EBHAHCED SORFACES 

Webb [Bef- 19], in an extensive review of the evolution 
of enhanced surface geometries, notes that the ability of 
roughness to improve nucleate boiling performance has been 
known for over 50 years. However, it was not until 1968 

that the first commercial enchanced surface was patented 
[Ref. 20]. Since then, the number of commercial enhanced 
surfaces has dramatically increased as the understanding of 
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their design and operation has grown. Of the many possible 
methods for heat- transfer enhancement, two areas are 
currently being commercially developed; 1) fins and surfaces 
with reentrant cavities, and 2) porous coatings. 

1 . Fins and Sur faces wi^ Ree ntrant C avitie s 

Ihe ability of surface abrasion, open grooves, and 
fins to improve the heat -transfer performance of a smooth 
surface was first studied in the 1930 *s. The main diffi- 
culty with using surface abrasion to improve the heat- 
transfer performance is that fouling of the surface 
eventually returns the performance to that of a non-abraded 
surface. Studies in the 1950’ s and 1960’s centered on fins. 
Recent comparative studies for refrigerants by Carnavos 
[Ref. 4] and Yilmaz and Hestwater [Ref. 2] found that fins 
and grooves result in a 50-100 percent permament improvement 
in the heat-transfer performance compared to a smooth plain 
tube in the same refrigerant. Hebb [Ref. 19], in his liter- 
ature review, describes how researchers in the early 1960*s 
found methods to improve the performance of fins by creating 
reentrant cavities on their surfaces. Reentrant cavities, 
such as shown in Figure 2.6, act as very stable nucleation 
sites and thereby enhance the heat-transfer performance. 
For a cavity to function as a nucleation site and remain 
active, even after the surface is subcooled, the mouth diam- 
eter (D) must fall within a critical range. Also, the 
cavity must have a reentrant shape with a maximum reentrant 
angle (6). The optimum mouth diameter (D) and reentrant 
angle (©) are functions of the fluid properties. 

The Gewa-T surface, patented in 1979 (manufactured 
by the Wieland Company), and the Ihermoexcel-E surface, 
patented in 1980 (manufactured by the Hitachi Company), are 
two commercial surfaces which use modified fin shapes to 
form the necessary reentrant cavities. Figure 2.7 shows the 
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(a) Schematic cross section of the Gewa-T surface 




(h) Schematic view of the Thermoexcel-E surface 



Figure 2.7 Surface Details of Gewa-T and 
Theraoexcel-E Reentrant Surfaces. 
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details of these tube surfaces. Carnavos [Ref. 4] found that 
in pure R-11, the Gewa-T surface outperformed a plain tube 
by 100-200 percent and the Thermoexcel-E surface outper- 
formed a plain tube by 300-400 percent. Curves 1-5 of 
Figure 2.8 show the relative improvement in the boiling 
heat-transfer performance of B-11 that can be achieved by 
using mechanically produced reentrant cavities. Work 
continues to optimiae these types of surfaces for the 
various refrigerants in use commercially. Both tubes have 
been tested in ref rigerant-oil mixtures and did not show a 
significant decrease in performance [Ref. 1 and 7]. The 
cost of these surfaces is not significantly higher than for 
smooth plain tubes, and the performance improvement is 

dramatic. 

2- P orous Coatings 

The second major type of enhanced surface is the 
porous boiling surface. Webb [Ref. 19] details the various 
production improvements and coating variations that have 
been made to the original 1968 patent by Milton of Union 
Carbide. The key to the performance of the porous coatings 

is their small reentrant cavities, which are interconnected 
by substrate tunnels. The particles used to make the coat- 
ings are usually copper or aluminum. According to Webb, 
researchers have found that the critical variable is the 
pore size rather than the particle size. Large pores are 
required for fluids with high surface tension and high 
thermal conductivity. Small pores are optimum for fluids 
with low surface tension and low thermal conductivity (like 
refrigerants). Curve 6 of Figure 2.8 shows the relative 
performance of the High Flux surface to finned tubes and 
mechanically produced reentrant surfaces. Carnavos [Ref. 4] 
found the High Flux surface to be 700-800 percent better 
than a smooth tube in R-11. No known studies have been 
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Figure 2.8 Comparison of R-11 Boiling Performance from 
Various Commercial Surfaces (from fief. 41 . 



published on the performance of porous coatings in 
refrigerant-oil mixtures. Some studies of the heat-tr ansfer 
performance of porous coatings in pure E-114 and 
refrigerant- oil mixtures have been made by Onion Carbide, 
but their results are not found in the open literature. 
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III. DESCHI PTION OF EXPERIMENTAL APPA RA TOS 
A. OVERALL APPARATOS 

An overall schematic of the experimental apparatus is 
shown in Figure 3.1, and a photograph is shown in Figure 
3.2. Karasabun [Ref- 8] describes the design, construction, 
and operation of the apparatus in detail. 

The apparatus consists of two Pyrex-glass tees. Liquid 
R- 1 1 4 is boiled in glass tee (1) and is condensed in glass 
tee (2) . Gravity drains the condensate from the condenser 
back to the boiling section. A water-ethylene-glycol 
mixture at -17 oc (1 °F) is pumped through the condenser 
cooling coil via a computer-controlled valve (VC) to 
condense the R-114 vapor. The sump (7) that supplies the 
water-ethylene-glycol mixture is cooled by a 1/2-Ton, E-12 

air-conditioning plant. 

Valve VC controls the R-114 liquid temperature and pres- 
sure. Figure 3.3 is a photograph of valve VC and the 
computer-controlled motor that operated VC. Opening VC 
causes more R-1 14 to condense and lowers the system pres- 
sure. Also, it returns more subcooled liquid to the boiling 
section which lowers the bulk liquid temperature. Since 
data at many heat fluxes was desired for a constant tempera- 
ture, it can be seen that changing the heat flux without 
adjusting VC would change the system pressure and tempera- 
ture. A computer-controlled valve was thought to be the 
best way to rapidly return the system to the desired satura- 
tion temperature following a heat flux change. Sections 
III.E and IV. D describe in more detail the computer- 
controlled valve and the operation of the system with the 
computer-controlled valve in use. 
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Schematic of Experimental Apparatus. 
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Figure 3,2 Photograph of Experimental Apparatus, 
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(b) Computer for Computer-Controlled Valve. 

Figure 3.3 Photographs of Computer-Controlled-Valve 

Eguipaent. 
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oil was added to the liquid R- 1 1 4 by draining it from a 
glass oil cylinder (5) . The oil cylinder was refilled as 
needed through valve 7-2 from the oil reservoir (4) . 

Two configurations of boiling tubes were tested. Short 
boiling tubes were tested to determine the correct assembly 
procedure to obtain data on the normal 431.8 mm (17 in.) 
long boiling tubes. The short tubes were cheaper to make, 
thus more debugging attempts could be made by testing them. 
Section III.C describes the details of the construction of 
the short and long test tubes. 



B. 



CIl SAHPLING APPAHATOS 



Following all data runs, an attempt was made to sample 
the local oil concentration in the vicinity of a boiling 
tube. Figure 3.4 shows the oil sampling apparatus. By 

opening valves S-1, S-2, and S-4, the probe line could be 
purged, trapping a sample inside a flexible silicon tube 
that was 30.5 mm (12 in.) long with a 3.16 mm (1/8 in.) 



3y weighing the sample 
leaving behind the oil. 



inside diameter using pinch clamps, 
tube and then boiling off the R-114 
the mass percent of oil in the R-1 14-oil mixture was deter- 
mined. The mass fraction of oil was calculated by: 
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A Pjrecisa Model 80 electronic mass balance was used to weigh 
the samples. The Precise Model 80 is accurate to 
± 0. 0001 g. 

C. BOILIHG TOBE CONSIBOCTION 
I • Short Tube s 

figure 3.5 shows the design of the short tubes. The 
short tubes were 15.9 mm (5/8 in.) in outer diameter, 
12.7 mm (1/2 in.) in inside diameter, and 203. 2 mm (8 in.) 
in length. The short tubes extended 152.4 mm (6 in.) into 
the liquid R-114 from the left end flange. A 25.4 mm (1 in.) 
long epoxy plug insulated the right end of the tube. A 1 mm 
(0.04 in.) thick copper disk behind the epoxy plug was soft 
soldered in place to act as a pressure barrier. The short 
tubes were heated by a 500-Batt 240-Volt stainless-steel 
cartridge heater. The heater was 6.35 mm (1/4 in.) in cuter 
diameter and 101.6 mm (4 in.) in length. 

The first short tube was made from thick-walled 
copper tubing. This tube was solid oxygen- free, high 
conductivity (OFHC) copper. Four 1-2 mm (3/64 in.) diameter 
holes were drilled into the wall of this tube at a diameter 
of 12.7 mm (1/2 in.) for thermocouple channels. Since this 
tube was solid, and had no sleeve interface, it did not have 
an interface resistance. Consequently, it was the reference 
tube against which all other tubes were compared to deter- 
mine the amount of contact resistance they had. 

Six other short tubes were made. Five were made of 
soft copper tubing and had sleeves inserted into the tube as 
indicated in Table 1. Soft soldering of the sleeves to the 
tubes was determined to yield negligible contact resistance 
by comparison with the solid tube. The last short tube (7) 
was made of 90-10 ccpper-nickel and was coated with High 
Flux over the active 101.6 mm (4 in.) long section. This 
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Figure 3.^1 Oil Sampling Apparatus. 




(a) Sectional view of short tube. 



A-A 




(b) Thermocouple sleeve unwrapped (at section A-A) to show 
the relative locations of the thermocouple channels 
(all dimensions in millimeters). 




(c) Left-end view of the short tube. 



Figure 3.5 Sectional ?ie«s of Short Boiling Tube. 
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tube tested the compatibility of the soft-solder-assembly 
method with the copper-nickel tube ' and High Flux coating. 
Section 7. A describes the results of data taken on the short 
tubes and the selection of the soft-solder method of 
assembly for the long tubes. 



TABLE 1 

Summary of Short Tube Assembly Methods 



Tube Surface 



Remarks 



1 

2 

3 

4 

5 

6 
7 



Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
High Flux 



Solid, thick-walled tube (reference) 
Slide-fit (0.005 in. clearance) 
Slide-fit (0.002 in. clearance) 
Press-fit (0.004 in. interference) 
Press-fit (0. 006 in. interf erence) 
Soft-soldered 
Sof t- solder ed 



2 . long Tube 

Figure 3.6 shows the design of the long boiling 
tubes. These boiling tubes were 15.9 mm (5/8 in.) in outer 
diameter, 12.7 mm (1/2 in.) in inside diamter and 431.8 mm 
(17 in.) in length. The center 203.2 mm (8 in.) was the 
active test section. For the copper-nickel tube, the center 
section was the only portion of the tube that was coated 
with High Flux. The remaining 114.3 mm (4.5 in.) on either 
side of the center section were smooth and unheated, and did 
not nucleate under any heat flux or oil condition. Karasabun 
[Ref. 8] describes hew these end-surfaces were treated by 
the data-reduction program as an extended fin from the 
center section and how their heat loss was accounted for. 

The center section was heated by a 1000-Watt 
240-Volt stainless-steel cartidge heater. The heater was 
6.35 mm (1/4 in.) in outer diameter and 203. 2 mm (8 in.) in 
length. The heater was surrounded by a copper sleeve with 
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Figure 3.6 Sectional View of Long Boiling Tube. 




eight 1.3 mm by 1.3 mm (0. 050 in. by 0. 050 in.) thermo- 
couple channels in them. Figure 3.7 shows the details of 
the channel layout. The thermocouple hot junctions were 
welded to the sleeve. Appendix A describes the calibration 
of the thermocouples . The channels were oriented to provide 
both axial and circumferential readings of the tube inner 
wall temperature. The sleeve was soft soldered to the tube. 
The data on the short soft-soldered tube closely matched the 
reference solid-tube data, and the short sof t-soldered-tube 
data matched similiarly with long tube data. The maximum 
circumferential wall temperature variation in the long 
smooth tubes was 0.80 K (0.31 °F) at 50 kW/m^ compared to a 
solid tube circumferential variation of 0.34 K (0-61 op) , 
Section V.A describes in more detail the circumferential 
variation of temperature that resulted using the various 
tube construction methods. Some axial temperature variation 
was experienced in the long tubes, particularly the long 
High Flux tube. Hon-uniform heat generation from the 
cartridge heater is believed to be responsible for the axial 
temperature variation. Section V.B descibes in more detail 
the long tube axial temperature variation. 

D. DATA ACQOISITIOH AHD BEDOCTION 

A Hewlett-Packard 3497A automatic data acquistion/ 
control unit was used to read thermocouple outputs and to 
read an analog signal representing the current and voltage 
supplied to the cartridge heater. A Hewlett-Packard 9826A 
computer unit was used to control the HP-34 97A and to 
analyze and store the data. 

Information was entered through the computer keyboard to 
prompt the HP-3497A to automatically scan each channel. All 
thermocouple measurements were accomplished by 0.245 mm 
diameter (30 gage) copper-constantan (type-T) thermocouples. 
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Top 

1 5 

I ^ 1 

4 X 2.6 X 8x 

I ^ 1 

Bottom 

(a) View of the boiling tube thermocouple locations as seen 
from the front of the experimental apparatus. 



A-A 




(b) Thermocouple sleeve unwrapped (at section A-A) to show 
the relative locations of the thermocouple channels 
(all dimensions in millimeters). 



LEFT-END VIEW RIGHT-END VIEW 




(c) End views of the boiling tube. 



Figure 3.7 Long Tube Thermocouple Channels 
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A power sensing device, described by Karasabun [Ref. 8],- 
converted the AC current and voltage values supplied to the 
cartridge heater to a 0-5 V analog, DC signal for scanning 
by the HP-3497A. Table 2 lists the channel allocations in 
the HP-3497A. 

Following data acquisition for each point, results were 
computed according to the step-wise procedure outlined by 
Karasabun [Ref. 8], and summarized in Appendix 3. Appendix 
B also includes a complete listing of the da ta-reduction 
program (DRP2) . 



TABLE 2 

HP-3497A Channel Allocations 



Channel 



Purpose 



33 

34 

35 

36 
62 
63 



Tube wall thermocouples (location T(1) 
to T (8) in long tube, T{5)-T(8) not 
used in short tubes) 

Liquid R-114 thermocouple T (9) 

BacXup liquid R-114 thermocouple T{10) 
R-114 vapor thermocouple T(11) 

Sump thermocouple T(12) 

Cartridge heater voltage analog signal 
Cartridge heater current analog signal 



E. CCMPOTER-COHTEOLLED 7AL7E 

The computer-controlled valve (VC) was a Whitney, 
screwed bonnet, regulating valve with 3/8 in. Swagelok 

fittings. The valve travel was 10.5 turns from full shut to 
full open. The valve handle was replaced by a 101.6 mm 
(4 in.) hard rubber disk (seen in Figure 3.3) which was 
rotated by the motor pinion gear. 

The computer-controlled motor was a General Electric 
"Minigear Motor” with a speed of 105.7 rpm and a torque of 
3.39 N-m (30 lb-in). The motor direction was controlled by 
two sets of Crydom solid-state relays that acted to open. 



shut, or hold the motor depending on signals from the 
computer. The computer controlled the amount of time that 
the valve was moving. Approximately 3 seconds were required 
for the motor to turn the valve one turn. 

The flow through the valve was checked by a flow meter 
and verified to be approximately linear (1 turn = 10 percent 
flow) for the flows most often used. The valve position was 
tracked by a 10- turn potentiometer, which was connected to 
the computer-controlled motor by another 101.6 mm (4 in.) 
hard rubber disk driven off the motor pinion gear. Hard 
rubber disks were used instead of metal gears to avoid 
damage to the gear teeth during the program debugging stage. 
The rubber disks allowed sufficient slipping when, for 
example, the computer sent a valve-open signal even though 
the valve was fully cpen. 

The E-114 liquid temperature input to the computer was 
provided by a separate copper-constantan thermocouple 
installed in the same liquid R-114 thermocouple well that 
the HP-3497A data acquistion/control unit used. The R-114 
ther omocouple emf for the computer was amplified by an Cmega 
thermocouple DC millivolt amplifier before input in to the 
computer. The data acquisition system and the computer- 
controlled valve system were completely independent systems. 

The computer used for the computer-controlled valve was 
an Octagon Systems SYS-2A microcomputer with an Esprit I 
terminal connected via an RS-232C serial port. Appendix C 
lists the control program used. The control program was 
written in NSC "Tiny EASIC. " The control algorithm simulates 
a proportional-integral-derivative (PID) controller to vary 
the valve opening and shutting times. In the final program, 
limits were placed on the numerical value of certain program 
variables to prevent register overflow, jamming the valve 
fully open/shut, and to lessen the impact of system noise on 
the response of the system. NSC Tiny BASIC is limited to 
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integers from -32000 to 32000 and fractional numbers are 
truncated. The values of the weighting factors (A,B,C) for 
the proportional term (E) , integral term (I) , and derivative 
term (D) were determined by trial-and-error. 

The control algorithm simplifies to the following lines: 



60 Input reauired temperature R 
230 Read E-1T4 temperature M 

260 Compute error (E), and change in error (D) 

262 Add 1 each loop to integral sum (I) if error is positive 
264 Add -1 each loop to integral sum fi) if error is negative 
400 Valve command V= (E/A) + (B*D) + (I/C) 

420 If (V>0) then open valve 
430 Otherwise shut valve 
540 GCIO 230 



Section V.C desrites in detail the system response under 
this algorithm. Computer control of valve VC was not used 
during the data taking as originally planned. Section V.C 
discusses the reasons for this decision. 
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IV. EIPEB IMEMTAL PBOCEDOfiES 

A. IHSTAllATION OP TUBE IH APPAEATOS 

Prior to installation in the apparatus, the boiling tube 
surface was cleaned with a 2 percent Nital solution (to 
remove surface oxidation and oil ) , rinsed with acetone, and 
air dryed. Short tubes tested with and without the above 
treatment showed no change in the heat-transfer performance. 
The treatment was effective in removing the slight surface 
oxidation present following soft soldering without changing 
the smooth or High Flux tube performance. Additionally, 
following testing with oil, a treatment was needed to return 
the High Flux surface to the ”no-oil" condition for further 
testing. 

After installing the tube in the glass tee, the appa- 
ratus was evacuated to 29 in. Hg by the portable mechanical 
vacuum pump (6) shown in Figure 3.1. System pressure was 
measured by a Marsh pressure gage (30 in. Hg to 150 psi 
range, ±2.5 in. Hg and ±0.5 psi accuracy). The apparatus 
was left at vacuum for 2 hours to check for leaks prior to 
each run. No noticable drop in vacuum was observed within 
the accuracy of the pressure gage. Next, the system gage 
pressure was raised to 0.19 MPa (27 psi), the saturation 
pressure of R- 1 1 4 at 21 oc (70 °F) , by opening valve V7 (see 
Figure 3.1 for the configuration of the valves). An 
Automatic Halogen Leak Detector, TIF 5000, was used to check 
for E-114 leakage. The sensitivity of this detector is 3 
ppm minimum concentration. After pressure egualization with 
the P.-114 reservoir (3), the reservoir drain valve (V6) and 
condenser return valve (V5) were opened to fill the boiling 
tee with liquid B-114. 
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Prior to installation, the left end flange had been 
scribe marked to indicate the licuid level corresponding to 
1600 cc (2500 gm) of liquid B-114 at 21 oc (70 op) . This 
was the mass of pure R-114 in the apparatus at the beginning 
of the data runs. The apparatus was now ready for taking 
data . 

B. G2NEEAL OPEEATIOH 

Table 3 lists the 109 data runs accomplished during this 
thesis effort and their purpose. The data runs were 
numbered sequentially and preceeded by a 2- or 3- letter 
prefix to indicate the tube type. The tube prefixes were: 

SS = Short Solid Tube 
HF = High Flux Long Tube 

WH = Eieland Hard Copper Long Smooth Tuba 

SSF = Short Slide- Fit Tube 

SPF = Short Press-Fit Tube 

SST = Short Soft-Soldered Tube 

SHF = Short High Flux Tube 

The short tube runs consisted of 6 data points at 6 
different heat fluxes (usually 59, 37, 22, 14, 8, and 5 

kB/m2) . The normal tube runs consisted of at least 7 
different heat fluxes with 6 data readings at each heat flux 
(usually 98, 6 1, 37, 22, 14, 8, and 5 kW/m2) with some addi- 
tional low heat fluxes investigated for the 0, 3, and 10 

percent oil cases to check for the onset of nucleate boiling 
and hysteresis in the High Flux surface. 

In all cases, the data set was begun by starting the 
cooling pump (8) and opening valve VC slightly to slowly 
cool the liquid R-114 to the desired saturation temperature. 
The R-114 vapor temperature was also monitored and when both 
liquid and vapor temperatures were stable (usually after 
30-45 minutes) , the heat flux would be established, the 
saturation temperature reestablished, and the data aquisi- 
tion unit was allowed to take data. Usually, the vapor 
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TABLE 3 

Summary of Data Runs 



Run 


Tsat 


Mo. 


No. 


(OC) 


Pts 


NH10 


SPF1 1 


6.7 


3 


SPF12 


-2.2 


6 


SPF13 


2.2 


6 


SPF14 


6.7 


6 


SSF15 


-2. 2 


6 


SSF16 


2.2 


6 


SSF17 


6.7 


6 


SSF18 


-2. 2 


6 


SSF19 


-2.2 


6 


SSF2 0 


-2.2 


6 


SSF2 1 


-2-2 


6 


SSF22 


6.7 


6 


SSF23 


-2.2 


6 


SPF24 


-2.2 


6 


SPF25 


-2. 2 


6 


SPF26 


-2.2 


6 


SPF27 


-2.2 


6 


SPF2 8 


6.7 


6 


SSF29 


-2. 2 


6 


SSF30 


6.7 


6 


SPF3 1 


-2.2 


6 


SPF3 2 


-2.2 


6 


SPF33 


2-2 


6 


SPF34 


6.7 


6 


SS35 


-2.2 


6 


SS36 


-2.2 


6 


SS37 


6.7 


6 


SPF38 


-2.2 


6 


SPF39 


-2.2 


6 


SPF40 


-2.2 


6 


SPF4 1 


-2.2 


6 


SS42 


-2.2 


6 


SST43 


-2.2 


6 


SST44 


-2.2 


6 


SHF45 


-2.2 


4 


SHF46 


-2.2 


5 


SHF47 


-2.2 


7 


SHF48 


-2.2 


8 


SHF4 9 


-2.2 


6 


SHF50 


-2.2 


6 


SHF 51 


-2.2 


6 


SHF52 


-2.2 


6 


SHF53 


-2.2 


6 


SHF 54 


-2.2 


6 


SHF55 


-2.2 


6 


SHF56 


-2.2 


6 


SHF57 


-2.2 


6 


SHF58 


-2.2 


6 


SHF59 


-2.2 


6 


SHF60 


-2.2 


6 


SHF6 1 


-2.2 


6 


SHF62 


-2.2 


6 


SHF63 


-2. 2 


6 


SHF64 


-2.2 


6 


SHF65 


-2.2 


6 


SHF66 


-2.2 


6 


SHF67 


-2.2 


7 


SHF6 8 


-2.2 


6 



Remarks 

Runs WH01 to WHIG taken by Karasabun 
Debug new tube 

Press-fit tube 4 (0.004 in. interference) 
Effect of Tsat 
Effect of Tsat 

Slide-fit tube 2 (0.005 in. clearance) 

Effect of Tsat 

Effect of Tsat 

Rotate tube 90 degrees 

Study low heat flux error band 

Study high heat flux error band 

Rotate heater -90 degrees, tube fixed 

Shift heater and thermocouples 

Repeatability 

Repeatability 

Repeatabili ty 

Shift heater and thermocouples 

Clean with Nital, acetone, and air dry 

Repeatability 

Slide-fit tube 3 (0.002 in. clearance) 

Effect of Tsat 

Repeatability 

Time Variation Study (1 day later) 

Effect of Tsat 
Effect of Tsat 
Short Solid Tube 
Repeatability 
Efiect of Tsat 
4 days later 

Shift heater and thermocouples 
7 days later 
10 days later 

Time Study Solid Tube (no effect) 

Short soft- soldered tube 6 
Repeatability 

Debug short High Flux tube 7 
Study Tvapor superheat 
Repeatabili ty 
Repeatability 

Study hydostatic head (+1 in. level) 
-Increase data pts at heat fluxes 
-Runs 49-54 form data set 
-Each run at different heat flux 
-Compare with set 55-60 
End of Set 49-54 

Study hydrostatic head (+0.5 in. level) 

-same as above set 49-54 

-same 

-same 

-same 

End of Set 55-60 

Install Tvapor radiation shield 

Study Tliquid subcooling 

Clean High Flux w/Nital and acetone 

Repeatability 

Repeatability 

Repeat ability 

Rotate tube +90 and +180 degrees 
Repeatability 
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Table 3 



Run 


Tsat 


No. 


No. 


(OC) 


Pts 


SHF 6 9 


-2.2 


7 


SHF70 


6.7 


36 


SHF71 


-2. 2 


36 


SHF7 2 


-2.2 


36 


SHF73 


-2.2 


36 


SHF74 


-2. 2 


36 


SHF75 


-2.2 


36 


SHF76 


-2.2 


36 


SHF77 


-2.2 


36 


NH78 


-2. 2 


9 


WH79 


-2.2 


42 


RH80 


-2.2 


1 1 


WH81 


-2.2 


42 


NH8 2 


-2. 2 


42 


NH83 


-2.2 


38 


HH84 


6.7 


49 


WH85 


-2.2 


66 


WH86 


-2. 2 


48 


WH87 


6.7 


48 


WH88 


-2. 2 


48 


WH89 


6.7 


48 


T?H90 


-2. 2 


36 


WH91 


-2.2 


48 


NH92 


6.7 


48 


WH93 


6.7 


48 


NH9 4 


-2.2 


48 


NH95 


-2.2 


61 


WH96 


6.7 


53 


NH97 


-2.2 


66 


1IH9 8 


-2. 2 


45 


NH9 9 


6.7 


66 


HF100 


-2.2 


66 


HF10 1 


-2.2 


72 


HF10 2 


-2.2 


24 


HF103 


6.7 


66 


HF104 


6.7 


69 


HF10 5 


6.7 


56 


HF106 


-2. 2 


48 


HF107 


6.7 


48 


HF108 


6.7 


48 


HF109 


-2.2 


48 


HFI 10 


6.7 


48 


HF1 1 1 


-2. 2 


39 


HF11 2 


-2. 2 


72 


HFI 13 


-2.2 


48 


HFI 14 


6.6 


48 


HF1 15 


-2.2 


48 


HF116 


-2. 2 


72 


HF1 17 


6.7 


72 


HF1 18 


6.7 


48 


SPF1 1 9 


-2.2 


36 


Note ; 


75 data ru 



2 percent 

3 percent 
6 percent 

10 percent 



Summary of Data Runs (cont'd) 



Remarks 

Apply thermal grease to heater (no effect) 

Effect of Tsat 

Repeatability 

Debug oil addition (0.2 percent added) 
Remove and clean tube, add 1 percent oil 
oil 
oil 
oil 

* oil 

Axial/Circumferential variation study 
0 percent oil, decreasing g 
Rotated tube -90 and -180 degrees 
Repeatability (of run WH79) 

0 percent oil, increasing a 
Repeatability (of run WH82) 

0 percent oil, increasing 
■ ■ ■ ’ decreasing 

decreasing 
decreasing 
decreasing 
decreasing 
increasing 
decreasing 
decreasing 
decreasin g 
decreasing 
increasing 
increasing 
decreasing 



0 

1 

1 

2 

2 

3 

3 

3 

6 

6 

10 

10 

10 



percent 

percent 

percent 

percent 

percent 

percent 

percent 

percent 

percent 

percent 

percent 

percent 



oil, 
oil, 
oil, 
oil, 
oil, 
oil, 
oil, 
oil, 
oil, 
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up to 2 K, 



temperature would read higher, up to 2 K, than the liquid 
temperature due to the vapor becoming superheated. The 
apparatus, though coded by the R-114, was still hot enough 
(60 op) to superheat the vapor. Measurements of the liquid 
temperature, system pressure, and vapor temperature (with 
vapor probe shielded by a radiation shield) confirmed this. 
The liquid temperature best indicated the saturation 
temperature. 

The HP-3497A data aguisition unit would scan each 
channel, compute the heat- transfer coefficient, print the 
results (an example printout is shown in Appendix D) , and 
store the data on the floppy disk. 

Following the taking of all data points, the data set 
was statistically analyzed by subroutine STATS to compute 
the average heat-transfer coefficient at each heat flux and 
the standard deviation of the 6 data points for a given heat 
flux. The standard deviation of the 6 data points was 
usually 0.5 percent for the heat flux and 1-2 percent for 
the heat-transfer coefficient. 

After all data sets at a given oil conceatra tion were 
complete, oil was added via valve VI. The oil immediately 
dissolved in the R-114. No carryover to the condenser was 
noted except for several small drops at 10 percent oil and 
the highest heat flux during the last few data sets. 
Foaming occurred with the addition of oil, and increased 
with both increasing heat flux and increasing oil concentra- 
tion. Figure 4. 1 to 4.3 show photographs at heat fluxes of 
30 kw/m2 and 98 kW/m^ for oil concentrations of 0, 3, and 10 
percent. When oil addition was not taking place, the oil 
cylinder and reservoir were isolated from the apparatus by 
valves V3 and V4 to minimize R-114 absorption by the oil. 
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(b) 1 percent oil at 100 kW/m^. Foaming begins to appear. 

Figure 4. 1 Photographs of Boiling and Foaming 
at 0 and i Percent Oil. 
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(a) 3 percent oil at 37 kW/m^. 




(b) 3 percent oil at 100 kW/m^. 

Figure 4,2 Photographs of Boiling and Foamino 

at 3 Percent Oil. ^ 
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(a) 10 percent oil at 37 kW/m^. 




(b) 10 percent oil at 100 kW/m^, 

Figure 4-3 Photographs of Boiling and Foaainq 

at 10 percent Oil. 
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C. SAflPlIHG OIL 



Oil sampling was done at room temperature to prevent the 
temperature of the sampling lines from boiling off the R-114 
and yielding false readings. The sampling procedure used 
was : 

1. Weigh the empty flexiile sample tubing and pinch 
clamps. Connect them to the sampling apparatus. 

2. Open S- 1 , S-4, and set S-2 to position 2 to provide a 
purge path for the sample flow. 

3. lower the sample container to a height below the 
level of the E-114 glass tee. 

4. After purging the sample probe and lines, use pinch 
clamps to isolate the 304.8 mm (12 in.) flexible 
sample tubing and trap a sample. Shut S-1 and switch 
S-2 to position 1. 

5. Disconnect the flexible sample tubing and immerse it 
in ice to lower its internal pressure below atmos- 
pheric pressure. 

6. Open S-3 to provide a vent for the purge container 
and pour the purge volume back into the boiling tee 
by lifting it to a height so that it flows by 
gravity . 

7. Weigh the sample line (the mass balance used was 
accurate to ± 0-0001 g) . 

8. Open a pinch clamp and allow air to warm the sample 
line and evaporate off the liquid R-114. Keep the 
sample line on the mass balance to collect any drops 
of oil that may splatter as the liquid R-114 evapo- 
rates. 

9. After allowing several hours for the R-114 to evapo- 
rate, reweigh the flexible sample tubing. 

10. Calculate the mass percent of oil using equation 
(3.1) . 
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V. SYSTEM OPE RATION AND PHOBIJM 



A. CIRCDHPERENTIAL WALL TEMPEBATOEE VARIATION 

Following the construction of the solid# thick-walled # 
reference tube# an investigation was begun to develop an 
assembly method which would result in negligible contact 
resistance between the inner tube wall and the copper 
sleeve# and yield data comparable to that for the solid 
tube. 

Short slide-fit tubes 2 and 3 (see Table 1) had sleeves 

that could be slid into the test tube with tube 3 having a 

tighter clearance. The slide-fit tubes exhibited a circum- 
ferential wall temperature variation of 7-14 K at 50 kW/m2 
in pure E-114 boiling at -2.2 oc. This variation matched 
data on a similar long tube (with slide-fit sleeve) tested 
by Karasabun [Ref. 8], Tube 2 had an average wall tempera- 
ture of 33 °C and tube 3 had an average wall temperature of 
25 oc at 50 kW/m2. The tighter clearance of tube 3 resulted 
in a lower contact resistance# and the wall temperature 

subseguently dropped. The short solid tube under similar 
conditions exhibited a 0.34 K circumferential wall tempera- 
ture variation and an average wall temperature of only 

10.7 oc, Sauer et al. [Eef- 21 J# with a similar experi- 
mental apparatus# used a mechanically-press-fitted brass 
sleeve for obtaining data. Tubes 4 and 5 had copper sleeves 
that were mechanically cold pressed into the tube. The 
diametral interference was 0.01 mm (0.0004 in.) for tube 4 
and 0.015 mm (0.0006 in.) for tube 5. The interface pres- 
sure obtained was calculated to be 15.2 MPa (2200 psi) for 
tube 4 and 22.8 MPa (3300 psi) for tube 5. During pressing# 
the sleeves were lubricated with glycerin ^ 

O O O 
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high-thermal-conduc tivit y compound that, according to 
Incropera and Dewitt [Bef- 22], should result in a 10-100 
times reduction in the contact resistance when combined with 
the contact pressures mentioned above. 

Tubes 4 and 5 had circumferential wall temperature vari- 
ations of 0.9-1. 5 K with an average wall temperature of 
16.2 oc for tube 4 and 14.7 °C for tube 5 during initial 
testing. However, with time, the circumferential wall 
temperature variation grew to 5-8 K and the average wall 
temperature increased to about 28 °C for both tubes. The 
resulting drop in heat-transfer performance is shown in 
Figure 5.1. This result is believed to be due to the phenom- 
enon referred to as ’’stress relaxation.” 

Stress relaxation is a form of creep. The Metal s 
H andboo k [Bef. 23] notes that copper alloys easily undergo a 
decrease in stress resulting from transformation of elastic 
strain into plastic strain in a constrained solid. The 
phenomena occurs even at relatively low operating tempera- 
tures (at 25 oc, an 80 percent drop in stress can occur in 
200 hours) , and is of most concern in applications like 
press-fits and solderless wrapped copper connectors. The 
time-dependent data of Figure 5-1 appears to be the result 
of a dropping press-fit interface pressure yielding a higher 
contact resistance, higher inner wall temperatures, and 
lower heat-transfer coefficients. 

Stephan and Hitrovic [Bef. 7] used a combination of 
mechanically press-fit and soft-soldered inner sleeves to 
obtain data on the Gewa-T surface in R-114. Soft-soldering 
was initially overlooked by this experimenter because the 
normal heating method used is an oxy-acetylene torch on the 
tube surface. The flame temperatures of an oxy-acetylene 
torch are above 800 oc (1500 OF) . The resulting oxidation 
and heat damage (the High Flux coating melts at 800 OC) to 
the High Flux surface was unacceptable. Tests showed 
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however, that the inner cartridge heater coaid be used to 
apply a controlled amcunt of heat to melt a low temperature 
(melting point 160 OC) solder, while closely monitoring the 
tube temperature using the sleeve thermocouples. This 
method was used to produce soft-soldered short tube 6 which 
showed circumferential wall temperature variations of 0. 8 K 
and an average wall temperature of 11.5 °C at 50 kW/m^. 
Figure 5.2 compares the results of the soft-soldered and 
solid tubes. The agreement between the tubes is excellent. 
The long boiling tubes matched the short soft-soldered data 
very closely. 

Tests made by rotating the test tubes 90 degrees and 180 
degrees showed that the slight circumferential wall tempera- 
ture variation of these tubes was due to the surface charac- 
teristics of the boiling tube rather than due to the 
thermocouples. Thermocouples located near more active 

nucleation sites of the smooth and High Flux surfaces had 
slightly lower local wall temperatures. The 0.8 K circum- 
ferential wall temperature variation that resulted from 
soft-soldering the sleeves of the tubes is much smaller than 
the 11 K variation reported by Sauer et al. [Sef. 21], and 
is about the same as the 1 K variation reported by Stephan 
and Hitrovic [Bef. 7], 

B. AXIAl TEHPEB2T0EE VARIATION OF LONG TDBES 

Karasabun [fief. 8] reported an axial temperature varia- 
tion of about 20 K along the inner wall of his slide-fitted 
sleeve. Stephan and Hitrovic [Ref. 7] reported an axial 
temperature variation of 1 K along the sleeve they used. 
The long soft-soldered tubes tested in this experiment 
exhibited an axial temperature distribution that varied with 
heat flux. Figures 5.3 and 5.4 show the axial temperature 
distribution as a function of position along the boiling 
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Figure 5.3 Axial Teaperatuie Variation of Smooth Tube. 
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surface for the smooth tube and High Flux tube respectively. 
The variation at 20 kW/m^ (the heat flux at which Stephan 
and Mitrovic cite their 1 K variation) is less than 1 K for 
both tubes. The maximum axial temperature variation is 3 K 
for the High Flux surface at 9 8 kW/m^. 

The cartridge heaters used were precision- wound , 
ceramic-core, magnesium-oxide-insulated, Incoloy-sheathed 
flATLOTi FIHEF.OD heaters. The heaters had 6.35 mm (0.25 in.) 
long lava-rock plugs at either end of the heater that 
reduced the actual heating length on each end. The heaters 
were initially believed to generate a uniform heat flux at 
all power settings, but the axial wall temperature data 
indicate the heat generation varied with the power level. 

Since the axial temperature distribution was fairly 
linear over most practical heat fluxes (less than 37 kW/m2) , 
the arithmetic average of all 8 wall ther mocouples was used 
to compute the heat- transfer coefficient. This results in 
slightly lower heat- transfer coefficients and is a conserva- 
tive estimate of the performance of the High Flux surface. 
Appendix E analyzes the resulting error in the heat-transfer 
coefficient from using the arithmetic average for calcu- 
lating the heat- transfer coefficient. 

C. CCMPDTEE-CONTROLIED VALYE 

Figure 5,5 shows the system response to closing valve VC 
manually and the improved response of the system when using 
the computer-controlled valve. 3y correctly cycling the 
valve between open and shut, the system saturation tempera- 
ture could rapidly be changed. Establishing a stable equi- 
librium temperature after a large valve movement was, 
however, difficult during the initial testing of the appa- 
ratus. The extensive trial-and-error testing for the 
correct weighting factors of the proportional, integral, and 
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derivative terms resulted in a large data base concerning 
the saturated system operating characteristics and a large 
amount of practice for the operator in manual control of 
valve VC following unsuccessful computer-ccntrolled 
transients. 

Manual operation of valve VC was done in conjunction 
with the section of the data- reduction program (lines 1890 
to 2485 of program DEf2 listed in Appendix B) that monitored 
the saturation temperature and the rate of change of satura- 
tion temperature continuously. It was found that, because 
of the better temperature resolution of the HP-3497A data 
acguisition unit, with manual operation and a trained oper- 
ator, temperature transients could be maintained within 
± 0.02 K. This was a much tighter control band than was 
possible with the computer-controlled valve because of the 
poor resolution of its temperature sensing input from the 
Omega thermocouple DC millivolt amplifier. Additionally, 
after much practice and experience, the manual method was 
found to be as fast as, or faster than the use of computer- 
controlled valve. 

To obtain the highest accuracy data in the shortest 
period of time, computer-control of valve VC was abandoned. 
Section VII. 3 includes recommendations for further improve- 
ments to the computer control system to restore its 
usefulness. 
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TIME <MCN5 

(a) Natural system response to shutting valve VC to cause 
a desired increase in system temperature. 




(b) Improvement in system response time using computer 
control of valve VC to reach a desired temperature 
(note the different time scale). 



Figure 5.5 System Response to a Transient with 
Com puter*“C on trolled Valve. 
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VI. BESPITS AND DIS CPSS ION 



A. BCILIBG PEBFOEMINCE OF SMOOTH TOBE 

Figure 6.1 shows the nucleate pool-boiling performance 
of the smooth copper tube in E-114. The behavior with no 
oil represents typical nucleate pool-boiling performance. 
From point A to point B, a continuous increase in vail 
superheat (T,, -T^^, ) is observed when heat flux is 
increased. No bubbles were observed in this region of the 

curve as this region represents natural convection. At 
point B, incipient nucleate boiling occurs. From point B to 
point C (or C for 10 percent oil), a reduction in wail 
superheat is observed while the heat flux is continuously 
increased. This region is known as the mixed-boiling 
region, where transition from natural convection to nucleate 
pool boiling takes place. In this region, the heated 
portion of the tube began to activate an increasing number 
of nucleation sites, while the unheated ends showed no 
bubbles. In fact, the unheated ends underwent only natural 
convection, due to axial conduction of heat along the tube 
wall, at all heat fluxes. The transition from natural 
convection to nucleate boiling occurred rapidly when there 
was no oil present. The surface would burst into nucleate 
boiling in less than a second after the first nucleation 
site became active. At point C (or C*), all the available 

nucleation sites were apparently active. After point C, the 
wall superheat again increases with increasing heat flux as 
shown in region C-D. In region C-D, no new nucleation sites 
were seen to become active. Instead, the bubble departure 
rate increased. When the heat flux is decreased after having 
established complete nucleate boiling, the curve follows a 
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Figure 6.1 Heat-Iransf er Performauce of Smooth Tube. 



different path (point D to point E) . The existence of 
stable nucleation sites, which remain active over a wide 
range of heat fluxes, results in better heat- transfer 
performance than in natural convection, resulting in a lower 
wall superheat. 

The effect of adding oil is, according to Stephan 
[Ref. 93 , to introduce a mass diffusion resistance and lower 
the heat-transfer coefficient. As seen in Figure 6.1, in 
region A-B, both the 3 and 10 percent curves were lower than 
the 0 percent oil curve. The 3 percent oil curve is lower 
than the 10 percent oil curve probably because of the non- 
linear physical property variations of refrigerant-oil 
mixtures. The non-linear variation of surface tension (see 
Figure 2.3) would not seem to be responsible for this 
anomally. The curves in region A-B support the contentions 
of C hongrungeong and Sauer [Ref. 10] and Thome [Bef. 14] 
that the non-linear variation of physical properties of 
refrigerant-oil mixtures, other than surface tension, 
explains the heat-transfer behavior of refrigerant-oil 
mixtures. The effect of adding 10 percent oil was to delay 
the transition to ccnplete nucleate boiling on the tube. 
With 10 percent oil, the surface developed patches of 
nucleation sites that spread slowly with increasing heat 
flux, until they covered the entire surface (point C'). 
Region E-E (or D'-E) shows that oil increased the wall 
superheat slightly for 3 percent oil and significantly for 
10 percent oil. Again, this agrees with the concept of an 
increased mass diffusion resistance by the addition of oil. 

Figure 6.2 shows the heat-transfer coefficient of the 
smooth tube in R-1 14-oil mixtures as a function of heat 
flux. The curves show the heat-transfer coefficient of the 
smooth tube in region D-E, after complete nucleate boiling 
has been initiated. The effect of adding less than 6 
percent oil is seen to be small (about a 10 percent 
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reduction in the heat-transf er coefficient). However, an 
oil concentration of 10 percent causes a more significant 
drop in the heat-transfer coefficient (from 0 to 35 percent 
depending on the heat flux) . 

Figure 6.3 shows more easily the degradation that oil 
causes in the boiling heat- transfer performance of the 
smooth tube. This figure plots the heat-transf er coeffi- 
cient of the smooth tube relative to the heat-transf er coef- 
ficient in pure R-114 as a function of oil concentration. 
The effect of oil can be seen to depend also on the heat 
flux. Oil can be seen to generally degrade the performance 
of the smooth tube, except at heat fluxes less than 5 
and oil concentrations between 2 and 8 percent. This 
behavior was also seen by Henrici and Hesse (see Figure 
2.2) . Since this curve shows the heat- transfer performance 
in region D-E, with complete nucleate boiling, tne non- 
linear variation of the physical properties of refrigerant- 
oil mixtures, including surface tension, again probably 
accounts for this anomalous behavior. Since no measurements 
of the physical properties of the R-1 14-oil mixture were 
made during this investigation, any possible non-linear 
property variations of the mixture used are unknown, 
requiring future work, 

B. BOILING PBRFORHASCE OF HIGH FLOX SORFACE 

Figure 6.4 shows the nucleate pool -boiling performance 
of the High Flux surface in R-1 14-oil mixtures. The small 
magnitude of the wall superheats obtained during nucleate 
boiling should especially be noted. The High Flux surface 
in pure R-114 showed typical natural-convection region (A-3) 
behavior. The incipient point (B) occurred at much lower 
heat fluxes and superheats than for the smooth tube as shown 
earlier. The transition to nucleate boiling (B-C) was 
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Figure 6.3 Relative Effect of Oil on Smooth Tube 
Boiling Beat-Transfer Performance, 
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similar to the smooth tube, occurring very rapidly (less 
than 1 second) . Nucleate boilng from the High Flux surface 
results in a drop of the wall superheat by about a factor of 
10. This is due to the extremely high density of nucleation 
sites present on the High Flux surface. From the low heat 
fluxes (1-10 JcW/m2) at which the transition to nucleate 
boiling occurs, it can be seen that the High Flux coating 
also assists in the activation of stable nucleation sites. 

Adding oil to the High Flux surface appears to delay the 
onset of nucleate boiling (point B wall superheat increases 
with increasing oil concentration) . However, the transition 
to nucleate boiling still occurred very rapidly, even at 1 0 
percent oil. The rapid transition to nucleate boiling (less 
than 1 second) on the High Flux surface is probably due to 
the interconnected cavities which can assist in nucleating 
the entire surface once a single site becomes initially 
active. Oil is unlikely to inhibit this characteristic of 
the High Flux surface though it apparently delays the 
initial activiation of the first nucleation site. As seen 
in region D-E, for heat fluxes less than 37 kW/m^, the 
effect of adding oil to the High Flux surface, once it has 
been nucleating fully, is to cause about a 30 percent 
increase in the wall superheat. At heat fluxes in excess of 
37 kW/m2 and at 10 percent oil, the wall superheat increased 
dramatically. As seen in Figure 6. 1, the wall superheat at 
a heat flux of 98 kW/m^ and 10 percent oil is about the same 
for both the smooth tube and High Flux surface. 

Figure 6.5 shows the heat-transfer coefficient of the 
High Flux surface as a function of heat flux for various oil 
concentrations. Again, oil is seen to degrade the nucleate 
pool-boiling heat-transfer performance. 

Air-conditioning plants typically operate in the heat 
flux range of 1 0 to 40 kW/m2 with less than 1 percent oil. 
In this region of pracical interest, the boiling 
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Figure 6.5 Boiling Heat-Transfer Coefficient 
of High Flux Surface. 



heat-transfer coefficient of the High Flux surface is about 
10 times better than a smooth tube. Also, it experiences 
only a 20 percent drop in the boiling heat-transfer coeffi- 
cient with the addition of oil (1 to 10 percent). This is 
seen clearly in Figure 6.6 which plots the heat-transfer 
coefficient, relative to the heat-transfer coefficient for 0 
percent oil, as a function of oil concentration. From 1 to 

10 percent oil at a heat flux of 14 kW/m2, the heat-transfer 
coefficient is about 80 percent of the no-oil heat-transfer 
coefficient. 

Figure 6.6 also shows that the oil-caused degradation of 
performance on the High Flux surface is nearly independent 
of oil concentration at practical heat fluxes. Only at a 
heat fluxof 98 kW/m^ and 6-10 percent oil, does the High 
Flux surface experience significant performance degradation. 
At high heat fluxes and oil concentrations, the effect of 

011 may be to ‘'clog” the interconnecting cavities of the 
High Flux surface due to boiling off of the R-114. Clogging 
the E-114 surface with oil would prevent replenishment of 
the nucleation sites with R-114 liguid, preventing the 
nucleation process and leading to higher superheats. The 
time-dependent behavior of the High Flux surface in high oil 
concentrations and at high heat fluxes was not studied in 
this experiment. 

C. CCHPABISON OF BIGH FLOX TO SMOOTH TUBE BOILING 

PEBFOBMANCE 

Figure 6.7 shows the heat-transfer performance of both 
the High Flux surface and the smooth tube as a function of 
heat flux. Again, the 7-10 times improvement in the heat- 
transfer coefficient by the High Flux surface is easily 
seen. At extremely high heat flux and high oil combinations, 
the High Flux surface performs only sligthly better than the 
smooth tube. 
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Figure 6.7 Coaparison of High Flux and Smooth Tube 
Boiling Beat-Transfer Performance. 



Figure 6.7 also shows that the heat-transfer curves of 
the High Flux surface and the smooth tube are not parallel. 
At low heat fluxes# the High Flux surface is more effective, 
in comparison to the smooth tube, in enhancing nucleation 
than at moderate heat fluxes. At high heat fluxes, the 
performance of the High Flux coating begins to converge 
toward the smooth-tube curve because the surface became 
vapor blanketed. This is expected since both tubes should 
perform about the same when they are completely vapor 
blanketed. 

Figure 6.8 shows the relative improvement of the High 
Flux surface over the smooth tube as a function of oil 
concentration. For the heat fluxes of practical interest in 
air-conditioning plants, the High Flux surface is 7-10 times 
better than the smooth tube. 

D. EFFECT OF SATOHATIOH TEMPEEATOEE ON BOILING PEBFOEHANCE 

As reported by Stephan [Eef. 9], the effect of 
increasing the saturation temperature for both the smooth 
tube and High Flux surface was increased heat- transfer 
performance. Figure 6.9 shows the improvement in heat- 
transfer performance in pure E-114 achieved by raising the 
saturation temperature from -2.2 OC (28 OF) to 6.7 oc 
(44 Of), high heat fluxes, little improvement is seen in 

the High Flux surface performance because the surface is 
nearly vapor blanketed with bubbles. 

Figure 6.10 shows the effect of saturation temperature 
on R-114-oil mixtures. Increased saturation temperature is 
again seen to improve the performance of the High Flux 
surface as well as the smooth tube, even with 10 percent 
oil. This is consistent with the no-oil results, tut 
contradicts Henrici and Hesse's data (see Figure 2.5). 
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Figure 6.8 Belative laprovement of High Flux 
Surface over Saooth Tube versus Oil. 
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E. 



lOCAl OIL SAflPLIBG EFFOfilS 



The results of the attempt to sample oil locally near 
the boiling tube were inconclusive. The method outlined in 
Section IV. C did not produce either repeatable or accurate 
results. Checks of the sampling method were made by 
sampling the bulk liquid with 0 and 10 percent oil at zero 
heat flux. The 0 percent oil check yielded oil concentra- 
tions from 0 to 2 percent. The 10 percent oil check yielded 
oil sample concentrations from 5 to 25 percent. 

Further refinements to the sampling apparatus are needed 
to permit an accurate local sample of oil in the vicinity of 
a boiling tube. Section 7II.B contains recommendations to 
improve the sampling apparatus. 
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VII. COMCI DSIONS AND RECOa HENDATION S 
COHClOSIOIS 

1. In pure E-114, the pool-boiling heat-transfer coeffi- 
cient of the High Fluz surface is about 10 times 
larger than that of a smooth tube. 

2. The High Flux surface began nucleate boiling at low 
heat fluxes, about 1 kW/m2, compared to a heat flux 
of about 10 kU/m2 for the smooth tube. 

3. Oil delayed the onset of nucleate boiling with the 
High Flux surface. 

4. Oil resulted in about a 20 percent reduction in the 
heat-transfer coefficient of the High Flux surface 
for most practical heat flux (less than 37 kW/m^) and 
oil combinations (less than 6 percent) . 

5. At heat fluxes of 98 kW/m^ and greater than 6 percent 
oil, the performance of the High Flux surface 
degraded by as much as 80 percent. The performance 
of the High Flux surface was little better than the 
smooth tube at 98 kW/m2 and 10 percent oil. 

6. The boiling heat-transfer coefficient of the High 

Flux surface was about 7 times better than that of a 

smooth tube for oil concentrations from 1 to 10 

percent over the range of heat fluxes employed in 
air-conditioning plants (10-40 kW/m2) . 

7. The boiling heat-transfer coefficient of both the 

High Flux surface and the smooth tube increased with 
increasing saturation temperature. The improvement 
decreased at high heat fluxes for the High Flux 
surface due to vapor blanketing of the surface. 



8. Local oil sample results were not accurate or repea- 
table enough to determine if a variation in the local 
oil concentration occurs near a boiling tube. 

B. EECOMHENDiTIOHS 

1. The present cartridge heaters should be replaced with 
more reliable heaters that will produce a uniform 
heat flux axially at all heat loads. 

2. A solid copper tube should be coated with High Elux 
and tested. Eue to the very small superheats of the 
High Flux surface, even a small amount of contact 
resistance could affect the data at low heat fluxes. 

3. The physical properties of the S-1 14-oil mixtures 
tested should be measured to obtain information to 
better explain the reasons for the heat- transfer 
performance of both the High Flux surface and the 
smooth tube. Particularly, the anomalous rise in 
heat-transfer performance of smooth tubes when 1 to 6 
percent oil is present in R-114 should be studied. 

4. A secondary heater should be installed in the boiling 
section to keep the total heat input constant. The 
heat input to the secondary heater should be 
increased when the heat input to the boiling tube is 
decreased and vice versa. This modification will 
maintain a constant heat load on the condenser and 
eliminate the need to operate valve VC, except to set 
the saturation temperature at the beginning of a run. 
Alternately, the effort to computer control valve VC 
should continue by obtaining an accurate thermistor 
with a high temperature resolution for use with the 
SYS-2A microcomputer. Improved temperature resolu- 
tion would allow the computer-controlled valve to 
operate properly and free the operator from the 
demands of manual control of valve VC. 
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The performance of a bundle of High Flux-coated tubes 
in E-114-oil mixtures should be studied. 

6. Data should be obtained on the High Flux surface over 
a wider range of temperatures and with oils of 
varying viscosities. 

7. The High Flux surface should be tested for time- 
dependent heat-transfer performance by taking data 
periodically over a long time during boiling. 

8. A small oil sample container with valves and a vent 
line should be manufactured. The flexible tubing and 
pinch clamps used in this experiment did not properly 
hold or vent the sample. 



83 



APPENDIX A 



THEBMOCOOPLE CALIBRATION 



Karasabun £Ref. 6] describes the thermocouple calibra- 
tion equipment in detail. Two thermocouples were calibrated. 
One was made from wire at the beginning of the roll, the 
other from the end of the roll, following the making of all 
thermocouples used in the apparatus and tests. 

Essentially, the manufacturer-supplied calibration equa- 
tion for the thermocouple wire, a seventh order polynomial, 
was corrected slightly by adding to it a second order curve 
fit of the variation of the manufacturer’s computed tempera- 
ture for a given emf from a known set of reference tempera- 
tures (measured using a Hewlett-Packard 2804A quartz 
thermometer with a temperature resolution of ± 0.0001 K and 
accuracy of ± 0.03 K) . 

The manufacturer's emf to temperature conversion equa- 
tion is: 



T = ag + a-]_E + a2E^ + agE^ + 

a e“^ + a E^ + a E® + a E^ 

4 5 6 7 

where 

T = temperature (OC) 
ag = 0.1 00860910 
a^ = 25727. 94369 
a2 = -767345.8295 
a3 = 78025595.81 
a^ = -9247486589 
ag = 6. 97688E+1 1 
ag = -2. 66192E+13 
a^ = 3. 94078E+14 

E = thermocouple reading (volts) 



(A. 1) 
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Figure A. 1 shows the guartz thermometer reading minus 
the thermocouple readings (discrepancy) versus temperature. 
The two thermocouples agreed to within 0.05 K of each other 
and the manufacturer's seventh-order polynomial needed about 
a 0. 1 K increase to more accurately convert the emf’s to the 
true temperature. The correcting second-order polynomial 
was: 



DCP = + b^T + b^TZ (A. 2) 

where 

DCP = discrepancy (K) 
b^ = 8.6268968E-2 
b^ = 3.7619902E-3 

b = -5. 0689259E-5 

2 

T = thermocouple reading 

(from eguation A. 1) (OC) 



Thus, the temperatures computed by the da ta-reduction 
program (DEP2) were emf's converted to temperature by egua- 
tion A. 1 with corrections for that temperature computed by 
eguation A. 2 added to the temperature to get the true 
temperature. 

Since the data-reduction program utilized differences 
between thermocouples in all computations, such as wall 
temperature minus saturation temperature, the corrections 
above were necessary only for the computation of items 
dependent on the absolute temperature, like the fluid prop- 
erties. Appendix E describes in detail the uncertainty 
analysis and the effect of wall temperature variation on the 
computation of the heat-transfer coefficient. Thermocouple 
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Figure A.l Thermocouple Discrepancy Correction. 



variations of the soft-soldered tubes were most 
to a slight amount of contact resistance and to 
characteristics of the tubes tested, since 
involving shifting the wall thermocouples did not 
data. 



likely due 
he surface 
data runs 
affect the 
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APPENDIX B 



DATA REDOCTION PEOGBAH 



The data reduction program below consists of the 
following sections: 



Main Program 
Sub Main 
Sub Plot 
Sub Poly 
Sub Flin 
Sub Stats 



Menu of subprogram options 

TaJce data or reprocess data 

Plot data on log-log scale 

Compute least-squares curve fit of data 

Plot data on ‘linear- linear scale 

Compute average and standard deviation of 



Subprogram Main consists of the following steps: 

1. Create data file for data and plot points. 

2. Select tube type. 

3. Monitor heat flux or saturation temperature to estab- 
lish steady-state conditions. 

4. Scan all channels listed in Table 2 and save in data 
file. 

5. Convert raw emf's to temperatures, current, and 
voltage. 

6. Compute the heat-transfer rate for the cartridge 
hater. 

7. Compute the average wall temperature, the wall super- 
heat, and the film temperature. 

8. Compute the physical properties of R-114 using given 
correlations at film temperature. 

9. Compute the natural-convection heat-transfer coeffi- 
cient of R-114 for the non-boiling ends of the test 
tube. 

10. Compute the heat loss from the non-boiling ends. 



da ta 
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11. Calculate the corrected heat flux from the test tube 
to the liquid B-114. 

12. Calculate the boiling heat-transf er coefficient of 
the E-114 from the test tube. 

13. Print data. Store heat flux and wall superheat values 
in plot file. 

The following is a listing of the complete data reduc- 
tion program (DRP2) written in Basic 3^0 for the 
Hewlett-Packard 9826 computer. 



I m r N/^Mt : imp 2 

1005* 0A1E; Octobe* 19, 1994 

1010* PE'.'ISEO: March 13 (985 

10151 

1020 fJELP 

1025 F’RINTEH I 

I0JO FRINf USING *4k '"Select ootion:""' 

1025 FRINT llSING "EX ’’0 ItUlng data cr r e-proc ea s i ng orevious data’"" 

1040 F'PINI USING ’ 6 X,**I IMotting iJala on Lcg^Log f 

104 5 PR INI USING "EX, “"2 Plot ting data on Linear"*" 

1050 PR INI USIriG '‘6X,"2 Mate cross- olot coefft file*’" 

1055 INFUl Irii) 

I0G0 IF ldo*0 IllL'N call Mein 

10G5 IF Ido-I IHEN CALL FIol 

1070 If Ido-J iHf U CALI PI in 

1075 IF Idu-3 iMEN CALL Ccef 

1 000 EMO 

1095 SUB ttjin 

1090 COM /Cc/ Ct 7 ) , leal 

1095 DIM E.-»f(l2) MI2) IHd(G) U2a( 6 ) .Oi aUi ) Doa< b > ,L a( fi ) ,Lua* 6 ) ,Kcua( R ) 

1100 UAIA 0. l0Ol)f>e9l . 25727 . 94269 - 767245 . 8295 . 7e02S>595 . 8 I 

1105 UAIm -924 74BG5B9 b . 9 /OdilC H I , - 2 . F 6 I 52F ♦ I 2 3.94 0 7BFM4 / 

1110 RFAU C( • > 

1115 |•IMNI^R IS 701 
1120 Cl FAR 709 
1125 PIER 

I I 10 IMFUI ’ENIFR MONIH 0/UE ANO 1 IME < MM ; PO : UH ; MM : S3 > ' ,0a t e» 

! I 35 0»nrUT 709» * IP" iLat. F 
1140 OUirui 70'4,-|O* 

I I 45 ENir.R 709;Ortlef 
1150 PRlnr 

1155 PR INI ■ Men 111 , ilate anrl lime :*i0riite^ 

1160 nmn 

1165 PR INI USING * 1 0 X , ' "Mil I r ; Irogr an none ; 0r\P2 

1170 PEEP 

1175 INFUl ’ENfli'R DISK NUMUER* Un 

1180 PRINT USING ' I 6 V . ■ *n i j.i nunbei - " * 27". On 

1105 PIER 

1190 iMPUr ENTfR INPUT MH»E ( 0- 1054 A , I =f 1 1 E > * , I n 
1195 Pirr 

1200 INPUT ’ENILR fMCRMOCOuPl.E TYPE ( 0 -Nl U , I *|>|. 0 ) ’ , I ca I I 

1205 IF In-'0 TUFll 
i:i0 PFlP 

1215 iNPUl "GIUE A NAME FUR IHf RAW DATA F U E * .02 ^ f i I eS 
1220 PRINT USING *l 6 X.’"New file nene: ’ " . I4A’ lOC^^f i I e< 

1225 PFFP 

1250 INFUl "INPIM SIZE Of I II E OUAl’.Si’el 
12 35 CRFhIE etiAf 02_filr5 S i rf I 
1240 ASSIGfJ 0file2 TO D2 Fi‘e» 

12151 

I250‘ OFMMt file until Nf iin FNf'UN 
1255 01., f I lei "DUMMY • 

1260 CREAIE BI'AI DI ,f 1 K-f :• i tI 
!2C5 ASSH N 01 I lei io ni f I le$ 

1270 OUlPLtT 01 I lei ; Dale'S 
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127S BEEP 

i:sa iMPur "giv/e a n^me for ime plot file" p_fiies 
1295 BEEP 

1Z90 INPUT "INPUT SIZE OF FILE BD/^r",5i3eZ 
1295 CREATE BOAT P._flle$ 5i:r2 
1J00 ASSIGN 0Plot TO P.nieJ 
1305 BEEF 

1310 INPUT "ENTER NUMBER OF OEFECTU'E TCS ( 0=OEFAULT > * , I dt c 
1315 IF ldtc“0 THEN 
1320 LdtcI-0 
132S Ldtc2«0 



1330 PRINT USING *16X,""No defective TCa e.'^lst"*'* 

1335 ENO IF 

1340 IF Idtc-1 THEN 

1345 BEEP 

1350 INPUT “ENTER OEFECTIUE TC LOCAT ION" ,L dt c I 

1355 PRINT USING "IGX/’TC defective at location "",0'iLdtcl 

1360 Ldtc2'0 

1365 ENO IF 

1370 IF Idtc-2 THEN 

1375 BEEP 

1390 INPUT "ENTER OEFECTIUE TC LOCATIONS" .Ldtc 1 ,Ldtc2 

13B5 PRINT USING ■1GX,""TC are defective at locations " " .0 ,4X ,0" sLdtc 1 ,Ldtc2 

1390 ENO IF 

1395 IF Idtc>2 THEN 

1400 BEEP 

1405 PRINTER IS 1 

1410 BEEP 

1415 PRINT "INUALIO ENTRY" 

1420 PRINTER IS 701 
1425 GOTO 1305 
1430 ENO IF 

1435 OUTPUT 0Fi lei iLdtcl ,Ldtc2 
I 440 • In* 1 00 I Ion 
1445 ELSE 
1450 BEEP 

1455 INPUT "6IUE THE NAME OF THE EXISTING OATA F ILE " ,02^f i 1 e* 

1460 PRINT USING ‘16X,""0ld file name: " " . 1 4A“ 1 02_f 1 I e$ 

1465 ASSIGN 0File2 TO 02^file« 

1470 ENTER $Flle2iNrun 
1475 ENTER §File2;0old$ 

1400 BEEP 

1405 INPUT "GU'E A NAME FOR PLOT F I LE " ,P _f 1 1 e< 

1490 BFEP 

1495 INPUT "INPUT SIZE OF FRE B0AT",Sire3 
1500 CREATE BOAT P^flleJ SizeZ 
1505 ASSIGN 0Plot TO P..file$ 

1510 PRINT USING ‘ 1 6X , * " Th l s data set ta^en on : * " , 1 4A" tOol d$ 

1515 ENTER @Fi le3iLdtc 1 .Ldlc2 
1520 IF Ldlcl 0 OR Ldtc2-0 THEN 

1525 PRINT USING " 1 6X . “ " Thernocouq le s were defective at locat Ions : * " ,2( 30 ,4X )“ i 
Ldtc I ,Ldtc2 



1530 

1535 

1540 

1545 

1550 

1555 

1560 

1555 

1570 

1575 

1590 

15B5 

1590 

1595 

1600 

1605 

1610 

1615 

1620 

1625 

1630 

1635 

1640 

1G4S 



EMO IF 

ENTER 9File2iltt 
ENO IF 
PRINTER IS 1 
IF rn*0 THfN 
DEEP 

PRINT USING MX , 

PRINT USING "6X, 

PRINT USING "6X . 

PRINT USING "6X . 

PRINT USING "6X . 

PRINT USING ‘•6X, 
f'RINT USING "6X , 

PRINT USING *6X . 

INPUT Itt 
IF 1 1 1 G THEN 
BEEP 

PRINT "INUALID ENTRY" 

GOTO 1560 
ENO IF 

OUTPUT aPtlel I Itt 
ENO IP 

PRINTER IS 701 

PRINT USING ■•|6X."-Tuhe Tvoe: 



"Select tube tyoe""" 

"0*Snocth 4 inch Ref""" 

"l*Snooth 4 inch Cu ( Fre s s / S I i do ) ‘ " * 
"Z*Soft Solder 4 inch Cu""" 

"3=>Soft Solder 4 u»ch HIGH FLUX’"" 
"4*Wieland Hard 0 inch*’" 

"5 UUGH PI UX B inch" " " 

"6»GEUA--K 19 F/in’“ 



0 '• I 1 1 1 
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I6S0 SEEP 

1GG5 INPUT 'ENrER OUTPUT UERSIDN ( 0=LDNG J -SHORT ,2*NDNE >", lov 
16601 

1GG5I Dl-Oirtneter at t hernocouo 1 e oositions 

IB70 DATA .0111 125 . .01 I I 126 01 I I 125, .0129540 ,. 01 2446 01 29540 ,.0100955 

IG75 READ DIaC' 

1GP0 01=D1a( Itt ) 

1685' 

1690' D2‘=Dlaneler of teat aection to the base of fina 

1695 DATA . 01 5875 01 5875 .. 01 5075 ,. 015824 0 1 5875 .. 0 1 5824 ,. 0 1 270 

1700 READ D2a( * ) 

1705 D2»02a(Itt> 

17101 

1715! Di*Inside dianeter of unefthanced ends 

1720 DATA .0127 . .0127 , .0127 , .0132 ,.0127 , .0132 . .01 I I 125 

1725 READ Dia( • ) 

1730 Di-Dia( Itt > 

1735! 

17401 Oo-Dutside diameter of unenhanced ends 

17 45 DATA . 0 1 5075 .. 0 1 5875 . 0 1 5875 ,. 0 1 5024 , . 01 5075 . . 01 5824 , . 0 1 270 
1750 READ Doa( • ) 

1755 Do»Doa(Itt) 

1760' 

1765' L=Length of enhanced surface I 

1770 DATA . 1016, .101E,. 1016, .1016, .2032, .2032, .2032 
1775 READ Lei( • ) 

1780 L=La<Itt> 

1785' 4 

17901 Lu=Length of uttenhanced surface at the ends 
1795 DATA . 0254 ,. 0254 ,. 0254 .. 0254 ,. 0762 .. 0762 .. 0 762 
1800 READ Lua( • ) 

1805 Lu-Lua<Itt> 

1 01 0 < 

1015! Kcu«Thernal Conductivity of tube 
1820 DATA 398 344.344,45,344,45.344 
1 825 READ Kcua( • ) 

1830 Kcu'KcuaUt t ) 

1835 A=PI»^Do 2-Di"2)/4 

1840 P=FI»Do 

1845 J-1 

1650 3x-0 

1855 5v-0 

1860 S<s=0 

1865 5..V-0 

1 870 Reoeat ; ' 

1075 ir ln=0 THEN 

1880 DN KET 0,15 RECOVER 1870 

1885 PRINTER IS I 

1090 PRINT USING "4X ."-SELECT OPTIDN--’ 

1895 PRINT USING "6X,""0-TAKE DATA-'*’ 

1900 PRINT USING •6X,--|-SET flEAT FLUX' •“ 

1905 PRINT USING “6X."''2 = SET Tsat*" 

1910 PRINT USING MV, ““NOTE: KEY 0 = ESCAPE*"" 

1915 BEEP 
1920 INPUT Ido 
1925 IF ldo=0 THEN 2495 
I 930' 

1935' LOOP TO SET HEAT FLUX 

1940 IF Irlo=l THEN 

1945 OUTPUT 709i"AR AF62 AL63 UR5" 

1950 PFEP 

1955 INPUT "ENTER DESIRED Qdo",0ado 

I960 PRINT USING ' 4X , " " DES I RED Qrlo ACTUAL 0do“"* 

1965 Ei-r=10O0 

197v3 FOR 1 = 1 TO 2 

1975 OUTPi.'T 709 » “AS SA“ 
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1 960 SuM-0 
1985 FOR Ji-I TO 5 
1990 ENTER 709iE 
1 995 Sun*Sun+E 
Z000 NEXT Jl 

2005 IP I-l THEN Uolt-$un*5 

2010 IF I»2 THEN Amp«E 

2015 NEXT I 

2020 Aado’VoM •Anp/(P1 ‘OZ^L) 

2025 IF ABS( Aqdo-Oqdo ) Err THEN 

2030 IF Apdo^Dodp THEN 

2035 BEEP 4000.. 2 

2040 BEEP 4000,. 2 

2045 BEEP 4000,. 2 

2050 ELSE 

2055 BEEP 250,. 2 

20G0 BEEP 250, .2 

2065 BEEP 250,. 2 

2070 ENO IF 

2075 PRINT USING * 4X .MZ . 30E ,2X MZ . 30E “ i Oqdp . Aqdp 

2080 WAIT 2 

2085 GOTO 1970 

2090 ELSE 

2095 BEEP 

2100 PRINT USING ” 4X ,flZ . 3PE , 2X .MZ . 30E " sOqdp . Aqdo 

2105 Err=500 

2110 UA I T 2 

21 15 GOTO 1970 

2120 END IF 

2125 ENO IF 

21301 

21351 LOOP TO SET T?al 
2140 IF Ido-'2 THEN 
2145 BEEP 

2150 INPUT “ENIER OESIREO Tsaf .Otld 

2155 PRINT USING OTaal ATfal Rate Tv Rate- 

2160 Oldl-0 
2165 OId2=0 

2170 OUTPUT 709»”AR AF33 AL35 UR5- 
2175 FOR I-l TO 3 
2180 Sun«0 

2185 OUTPUT 7091-A5 SA’ 

2190 FOR Jl=l TO 2 
2195 ENTER 709.EIiq 
2200 Sun*5un+Ello 
2205 NEXT Jl 
2210 EliQ=Sun/2 
2215 Tld'FNTvav( FI in) 

2220 IF r=l THEM Atld»Tld 
2225 IF 1-3 THFN Tv=Tld 
2230 NEXT I 

2235 IF ABSf At Id-Otld )'.2 THFN 

2240 IF At Id'Ot Id THEN 

2245 PEEP 4000,. 2 

2250 BEEP 4000,. 2 

2255 PEEP 4000,. 2 

2260 ELSE 

2265 BEEP 250 . .2 

2270 BEEP 250,. 2 

2275 PEEP 250 .2 

2280 end ir 

2295 Err I =* At I d-0 1 dl 
2290 Oldl'Atld 
2295 Efr2*rv-Old2 
2300 Old2''Tv 

2305 PRINT USTNG " 4 X . 5< MOO . 00 I x ) - , 0 1 1 d A t I d Err 1 . T v . Err 2 
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2310 WAIT 2 
2315 GOTO 2170 
2320 ELSE 

2325 IF. A8S(Atld Ot ld)>. 1 THEN 

2330 IF Alld’Otld THEN 

2335 BEEP 3000 .. 2 

2340 BEEP 3000,. 2 

2345 ELSE 

2350 BEEP 800,. 2 

2355 BEEP 800,. 2 

23B0 END IF 

2365 Errl-Atld-Oldl 

2370 Oldl-Atld 

2375 Err2=Tv-Dld2 

2380 Old2»Tv 

2385 PRINT USING “ 4X ,5( NOD. DO , 1 X )" i D 1 1 d , At Id , Err I , Tv ,Err2 

2390 WAIT 2 

2395 GOTO 2170 

2400 ELSE 

2405 BEEP 

2410 Errl -At Id-Dldl 
2415 Oldl-^tld 
2420 Err2=Tv-Dld2 
2425 Dld2»rv 

2430 PRINT USING “ 4X ,5( NDO . DO . I X ) * i D t I d , A t 1 d .Err 1 , Tv ,Err2 

2435 UAIT 2 

2440 GOTO 2170 

2M5 END IF 

2450 END IF 

2455 END IF 

2460* ERROR TRAP FDR Ido OUT OF BDUNOS 

24B5 IF ldo'‘2 THEN 

2470 PEEP 

2475 GOTO >890 

2480 ENO IF 

24851 

2490! TAKE DATA IF lm-0 LOOP 
2495 PFEP 

2500 INPUT -ENTER BULK OIL Boo 
2505 OUTPUT 709i"AR AF25 AL36 URS* 

2510 FOR i-i ID 12 
2515 OUTPUT 709 r* AS SA" 

2520 Sun-0 
2525 FDR Ji-1 TO 2 
2530 ENTER 709«E 
2535 5iifn = S(jfn + E 
2540 NEXT li 
2545 Emf( I )-Sun/2 
2550 next I 

2555 OUTPUT 709i"AR AF62 ALB3 UP5* 

2560 FOR I-l TO 2 
2565 OUTPUT 709; “AS SA“ 

2570 5un=0 

2575 FOR Ji-1 TO 2 
2580 ENTER 709 lE 
2585 9uri-Sufn.tE 
2590 NEXT Ji 

2595 IF l»l THEM Ur-Sun/2 
2600 IF 1-2 "HEN Ir-Sun/2 
2605 NEXT I 
2610 ELSE 

2615 ENTER PF 1 1«2 ; Boo , To 1 dS ,Enf ( ♦ > , Ur , Ir 

2620 ENO IF 

2G25! 

2G30I COM'.'ERT enf'S TO TEMP , HOLT .CURRENT 
2635 Tmc»-0 



I 
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2G40 


FOR I-l TO 12 




ZG45 


IF ldtc;0 THEN 




2S50 


IF I-Ldtcl OR I“Ldtc2 THEN 




ZG55 


T< I )=-99.99 




ZSGO 


GOTO 2710 




2GG5 


ENO IF 




2S70 


END IF 




2675 


IF Itt«4 THEN 




2600 


IF I>4 mo I<9 THEN 




2605 


T( I )--99.99 




2690 


GOTO 2710 




2695 


ENO IF 




2700 


ENO IF 




2705 


T( I )»FNTvsv(Enfn 1 ) 




2710 


NEXT I 




2715 


IF Itt<4 THEN 




2720 


FOR I-l TO 4 




2725 


IF I=Ldtcl OR I“Ldtc2 THEN 




2730 


T Uid*T Wd 




2775 


ELSE 




2740 


T Wd^Tud^ T ( I ) 




2745 


ENO IF 




2750 


NEXT I 




2755 


Twd=Twd/( 4-Idtc ) 




2760 


ELSE 




2765 


FOR I-l TO 9 




2770 


IF I-Ldtcl OR I'Ldtc2 THEN 




2775 


Tuid*Twd 




2780 


El SE 




2785 


T tJd^ f ujd ♦ T ( I ) 




2790 


ENO IF 




2795 


NEXT I 




2 000 


Tud *Tud/( 8 I dtc .) 




2005 


END IF 




2010 


Tld’T( 9) 




2015 


Tld2=T( 10) 




2020 


Tv-T< 1 1 > 




2925 


IF Ht 3 THEN 




2830 


Tld2--99.99 




2035 


Tv-(T(l0)+T(ll)>/2 




2040 


END IF 




2845 


T Slino ’ f f 12) 




2950 


/>niJ * I r 




2855 


Uclt-Hr.:5 




2860 


0 = 1^01 1 'ftmo 




2865 


ir ltt*0 THEN 




2870 


Fc'j=FMKcu< T wd ) 




2875 


ELSE 




2000 


FcU'Kcua ( H t ) 




2885 


ENO IF 




2890' 

2095' 


FOURIER CONO'JCTION EDUOTION UlTH 


CONTACT RESISTANCE NE6LECTE0 


2900 


lu = Tud-0*L0G<02'01 )/(2*PI »Kcu»L ) 




2905 


Thctab-Tw-T Id 




2910 


IE Thetdb 0 THEN 




2915 


GEEP 




2920 


INPUT "TUALL T5AT ( O-rnuT INl.'E , 1 


“END ) " , I e V 


2925 


IE lev»0 THEN GOTO 1875 




29ZO 


IF I-v-1 TIIFN 3330 




2935 


ENO IF 




23401 

2945' 


COriPUIE ^'ARIOl.'S PROEFRTIE? 




2950 


Tf iln^ENFft lm( TuTId) 




2955 


Rho*FNPho ( T f i 1 fn ) 




2 360 


Mu = EMMu( T f 1 In ) 




29G5 


K“FNK( Tf 1 U ) 





9.4 



2970 Co-FNCo( Tf iln) 

2975 Beta-FNBet 0 ( Tf 1 Im ) 

2980 Ni“Mu/Rho 
2935 Aloha-K/(Pho»Co ) 

2990 Pr«Ni/Alpha 

2995 Psa^'FNPsaK Tld) 

3000' 

3005' COMPUTE NATURAL -CONUECTIUE HEAT-TRANSFER COEFFICIENT 
3010' FOR UNENHANCEO ENO( S ) 

3015 Mbar=l90 

3020 Fft-(Hbar»P/(Kcu»A) .5«Lu 
3025 Tanh-FNTanh(Fe ) 

3030 Theta-Thetab»Tanh/Fe 

3035 X •.“( 9 . 8 1 «Oe I a • The I ab *00 3 • Tanh/ ( Fe*N i • A loha ) T . I GGBS7 
304 0 Yy = ( I+( .559/Prl '(9/lB ' ) ( 0/27 ) 

3045 Hbarc-K/Do*( . B + . 3B7 • X... /Yy )2 

3050 IF ABS( (Hbar-Hbarc )/Mbarc .001 THEN 

3055 Hbar*( Hbar*Hbarc )* .5 

3OG0 GOTO 3020 

3055 END IF 

Z070! 

3075' COMPUTE HEAT LOSS RATE THROUGH UNENHANCEO END( S ) 

3080 01»(Mbar*P*Kcu»A)'.5*Thetab*Tanh 
3085 0c-0-2*0I 

3090 As=PI«D2*L 
3095' 

3100' COMPUTE ACTUAL HEAl FLUX AND SOILING COEFFICIENT 
3105 Odo'Oc/As 
3110 Htubc*Odo/ That ab 
31 151 

3120' RECORD TIME OF DATA TAKING 
3125 IF ln=0 THEN 
3130 OUTPUT 709i“TD“ 

3135 ENTER 709 i Told* 

3140 END IF 
3145' 

3150' OUTPUT DATA TO PRINTER 
3155 PRTNTFR IS 701 
3 1 BO If lov-0 IIIFN 
3IG5 PRINT 

3170 PRINT U5IMG -10X,"“Oata Set Munber - “ * .ODD . 2X . ‘ ‘Bu 1 1 Oil X - “ “ , DO. 0 ,5X , I 
4A ' I J ,Bod , T o I d$ 

3175 PRINT 

3180 PRINT USING " 1 0X , * " TC No : I 2 3 4 5 G 7 



3195 PRINT USING " 1 OX , “ Ter.o : “ " . 6 ( 1 X .MOD . DO ' * i T ( I ) . T( 2 > , T( 3 ) . 1 ( 4 ) , T ( 5 ) . T ( G ) . T ^ 
7 ' ,T( 8 ) 

3190 PRINT USING " 1 0X . " Tui.i lllod TliqdC T.’aor Psat Tsunp"**" 

3195 PRINT USING " 1 OX ,2 ( MHO . DO . I X ) , I X ,flOD . DO . I X . 2 ( I X .MOD . DO ) ,ZX ^ MOO . O' 5 Twa , T I d , 

T Id2 ,Tv ,Psa I ,Tsuno 

3200 PRINT USING " 1 0X . ' ” Thetab Htube Odo"“" 

3205 PRINT USING “ 1 0X ,MPD . 30 . I x . m2 . 30E , 1 X ,M2 . JOE *' i The t ab . Htube ,0do 

3210 END IF 

3215 IF Iov=l THEN 

3220 IF J=l THEN 

3225 PRINT 

3230 PRINT USING "10X.“'*RUN No OiIX Tsat Htube Odo Thetab‘"“ 

3235 END IF 

3240 PRINT USING '12X.30 4X DO .2X .MOO. DO , 3( I X ,M2 . JOE ) “ i J ,8op . T1 d .Htube ,Qdo .The t 
ab 

3245 END IF 
3250 IF In-0 THEN 
3255 BEEP 

32G0 INFI.'I 'OK TO STORE THIS DATA SET < I »Y 0-N ) 7 •• .01 
32G5 END IF 

3270 IF 01''! OR Im^I THEM J*1H 



I 



as 



3275 IF 0K»! AND ln=0 THEN OUTPUT @Fi le I iBoo . TcIdS Enf< • > , Ir 

3380 IF In=l OR 01=1 THEN 0U1FUT 9PlotiOdo Thetab 
3285 IK ln=0 THEN 
3230 BEEP 

3395 INPUT ‘WILL THERE BE ANOTHER RUN <1 - Y ,0-N )? ” .Go^on 
3300 Nrun=J 

3305 IF Go_onv'l THEN 3330 
3310 IF Go_on*l THEN Reoeat 
3315 ELSE 

3330 IF J-Nrun+1 THEN Reoaat 
33C5 ENO IF 
3330 IF In-0 THEN 
3335 BEEP 

3340 PRINT USING " 1 0X , “ “fJOTF ; "",ZZ,““ data runs were stored in file ■*,I0A*|J- 
I ,D2_f i J e? 

334«S ASSIGN @Filel TO • 

3350 OUTPUT @File2«Nrun-l 
3355 ASSIGN SFilel TO Ol.filet 
3360 ENTER @F i 1 e 1 . Oa t ef .L d t c Mdt c3 . It t 
3365 OUirUI (»Fi I rTiDatei L d I c I .Ldt r 2 . 1 1 1 
3370 FOR 1=1 TO Nrun-I 

3375 ENTER t?F i I e 1 i Boo . To 1 d« En f ( • ) .Ur . Ir 
33R0 OtiTPUT @Fi le2«Boo ,To) di CnfO.Ur.Ir 
53P5 NEYT I 

3390 ASSIGN i«F i lei TO • 

3395 PURGE “OUMNY* 

3400 ENO IF 
340S BEEP 

3410 PRINT USING 'l0X,--NOrE: “*.22.“ Y y oairs were stored in plot data file 
10A“iJ-l P^PlleS 
!4 IS AGSIGM t»Ei le2 TO • 

3420 ASSIGN 0riot 10 • 

3425 CALL Stats 
3430 8FEP 

3435 INPUT -| IHE TO PLOT OATA ( I = Y ,0-N 1 ? " ,0t 

3440 IF 01 »1 THEN 

3445 CALL Plot 

34SC END IF 

3455 SUBKNfJ 

3460 I 

34651 CUR‘'E KlIS OF PRUPFRTY FUNCTIONS 
3470 OFF FNKcut I ) 

3475' OFHC COPPER 250 TO 300 K 

3480 T1 =T+273. 15 'C TO 1 

3485 K-434- . 1 12»U 

3490 RETURN K 

3495 FNENO 

3500 OFF FNNu(T) 

3505' 170 TO 360 K CURVE FIT OK VISCOUS Ilf 

3510 T1 «1 ^273. IS !C TO K 

3515 Mu“EXP( -4 . 4636*( 101 1 . 47/Tl ))*1.0E-3 

3520 RETURN Mu 

3525 FNFMO 

3530 OFF FnCo(T) 

35351 IPO TO 400 K CURVE FIT OF Co 
3540 n = T »273. IS 'C TO K 

3545 Cp- . 4018BH .G5007E-3»ri f I . 5U94F-6*TT ~ 2-6 . G7R53E- 1 0» Tk 3 

3550 Co-Co ‘1000 

3555 return Co 

35G0 FNFNn 

3565 DfF FtlRt,o(T) 

3570 Tl •*r073. 15 1C TO 1- 
3575 A = 1 -M .8*n /7S3.9S) 'K TQR 

35P0 Ro = 36.32'6l . I46414 M ( l 3 > H 6 . 4 1 80 1 5 M H 7 . 4.76 838 • K " . 5 H . 1 ig028«K '2 
*585 Ro»Ro« .062429 
3590 RETURN Ro 



I 



as 



3595 FNEND 
3GP0 DEF FfJPr( T > 

3G05 Pr»FNCp< T )«FNMu( T )/FNK( n 
3B10 PFIURN Pr 
3615 FNEND 
3630 DEF FNK(T) 

36351 T:360 K WITH T IN C 
3G30 K*. 07I-. 000361 *T. 

3635 RETURN K 

3640 FNEND 

3645 DEF FNranh(X) 

3E50 P“EXP(X) 

3655 0-l/P 

3660 Tanh=<P“0 >/(P+0 ) 

3665 RETURN Tanh 

3670 FNEND 

3675 OFF FNT-'S'.(U) 

36P0 CDI-1 aCc/ C{7),Ical 
3585 T=HO) 

3590 FDR 1=1 TD 7 
3635 T = TtCU)*U-I 
3700 NEXT I 
3705 IF Ical=l THEN 

3710 T-T-6,74?3g34E-2fT*(g.0377043E-3-T*( 9.33539I7E-5) ) 

3715 ELSE 

3730 T = T ^8. G368g7E-2 + T*< 3. 76 tg'3E-3- 7*5.06893595-5) 

3735 END IF 
3730 RETURN T 
3735 FNEND 
3740 DEF FNBela( T ) 

3715 Rco=FNRho< T+ . 1 ) 

3750 Ron-FNPhoI 7-. 1 ' 

3755 8eta»-37(Roo + Ron)*(Roo-Ron )/.3 
3760 RETURN 8s t a 
3765 FNEND 

3770 DEF FHTf iln( Tw Jld) 

3775 Tf 1 ln“( Tw^ Tld'^3 
3780 RETURN Tfilm 
37P5 FNFN(1 
3790 OFF FMRsatfTc) 

3795! 0 TO 80 dsg F CURUE FIT OF Psat 
.3600 Tf=|.8*Tc + 33 



3805 

3810 

3815 

3G30 

3835 

3F30 

3835 

3840 

3645 

3?'-^0 

3855 

3P60 

3865 

3870 

3875 

3080 

3885 

3090 

3895 

3900 

3905 

3910 

3915 

3930 



Pa-5. 945535 + Tf*( . 1 5353083 f T f • ( 1 . 4840963E-3+ T f * 9 . G 1 5067 1 E-G ) ) 
Pg-Pa-1 4.7 

IF Pg 0 THEN ' +»P5I6,- = in llq 

Psdt'Pg 

FLEE 

P5at=Pg*29. 93 <• 1 4 , 7 
END If 
PETURN Psat 
FNEND 
SUB Plot 

CDM /CoW/ AU0.10) CUO) B( 4 ) ,Nop , Iprnt ,0po , I log 
INTEGER Ii 
PRINTER IS I 
8FrP 

INPUT "1 IKE DEFAIJl T HAI UFS FOP PI DT ( 1 - f , 0-N )? ” , Idv 
BEEP 

PRINT USING ”4x . "-S-l »*rl! notion:"’* 

FRTNT USING “6X **0 p versus delta-T"** 

PRINT USING ••GX,’*1 h versus delta-T*"" 

PRINT USING "RX,‘"3 h versus n"** 

INPUT Onn 
PFfP 

INPUT "SELFCT UNITS < 0=S I , 1 NGL ISM ) " Inn 
PRINTER IS 705 



3^:5 IF Idv;'l THEN 

BEEP 

INPUT "ENTER NUMBER OF CYCI ES FOR V-A!<I5’,C.< 
3340 BFEP 

3343 INPUT “ENTER NUMBER OE CYTLEB EOR Y-AXIS*, Cy 
3350 PEEP 

3353 INPUT 'ENTER MlN X-UALUE ^MULTIPLE OE 10)* Xntn 
3360 PEEP 

3965 INPUT "ENTER MIN Y-UALUE (MULTIPLE OF l0)*.Ynin 

3370 ELSE 

33?5 IF Ooo=0 THEN 

33R0 Cy=3 

39S5 C.x=3 

33RO Xp»in=.1 

3333 YniM»IO0 

4000 END IF 

4003 IF Ooo=l THEN 

40)0 Cy*3 

4015 C.^-3 

4030 Xmin^.l 

40ZI5 Yni n= I 00 

40^0 END IF 

aO'^S IF npc*2 THEN 

4040 Cv=3 

4045 c-=: 

4050 Xm I n* I OOO 
4 055 Yn i 100 
40F0 ENO IE 
4065 Efin (E 
4070 BEEP 

4075 PRINT “iNtSPliIP 7300 ::O0 07OO,SflO0i’ 

40P0 PRINT *SC 0,100 0.)00»IL 2.0:* 

40P5 Sf.-100/C,. 

4050 Cfv-IOO/Cv 
4035 PRINT "PIJ 0,0 PO* 

4)00 Nn=3 

4)05 FDR 1=) TO C, M 
4J 10 Y^r-<nin* )0' ( I- ) ) 

4) )5 IF 1=0. H THEN Nn») 

4)20 FOR .) = ! TO f)n 

4)75 IF J=) niEN PRINT " T| 2 0“ 

4)30 IF .1 = 2 THEN PRINT “Tl ) 0* 

4)35 Xa-<.9t*) 

4)40 XaiOT(Xa/Xmin)*5f< 

4 145 PRINT *PA" ;X . " 0i XT j * 

4)50 NEXT J 
4)55 NEXT I 

4)60 PRINT *PA ) 00 0:PUi " 

4)65 PRINT *PIJ PA 0.0 PO" 

4)70 Nn-3 

4)75 FOR 1=) TO Cy>) 

4)P0 Ya t = Yn : n • ) 0 " ( I - I ) 

4) PS IF I=CvH T)IFM Nm-I 

4)50 FOR J* ) TO Nn 

4)35 IF ) = ) T))EN f’RINT “ T( 2 0 

4200 IF .1 = 2 IHFN PRINT *TI ) 0’ 

4203 r>^»Yat*I 

42)0 Y = l FT{ Ya/Ynjn ''Sfy 

42)5 PRINT *PA 0 * I Y *VT“ 

4220 NEXT 
4225 next I 

4730 PRINT "PA 0 )00 TI. 0 .3* 

4715 Nn=3 

4240 FOR 1=) TO C. M 

4245 X,,t = Xn in» ) 0 < I - ) » 

4 250 IF I 4 ) THFN Nn = ) 



as 



4:?55 FOR J-1 TO Nn 

42G0 IF J=I THFN PRINT “ Tl 0 2” 

42BS IK J'i I THEN PRINT "TL 0 I “ 

4270 = 

427S X-1 GTCXa/Xnln )*Sf^ 

4200 PRINT "PA“ iX , , 1001 XT" 

42fiS NEXT J 

4290 NFXT I 

429S PRINT 'PA 100,100 PU PA 100 0 PO" 

4300 Nn*9 

4305 FOR I-l 10 Cy+1 
4310 Yat"Ynin*IO''(I-l) 

4315 IF I-Cy+I THEN Nn- 1 

4320 FOR J-1 TO Nn 

4325 IF J=1 THFN PRINT " TL 0 2' 

4330 IF J>1 THFN PRINT 'TL 0 I" 

4335 Ya-Yat*J 

4340 Y-LGT( Ya/Vnin >*Sfy 

4345 PRINT “PD FA 100,' Y,”YT“ 

4350 NEXT J 
4355 NFXT I 

4360 PRINT 'PA I P0 . I 00 PIT 
43B5 PRINI "PA 0 -2 ER 1 .5 2’ 

4370 Ii-lGT(Xnin) 

4375 FOP I-l TO n.*l 
4390 Xa--Xnin* 10 'f I-l > 

4395 X=LGT(X^/Xntn >*5f> 

4390 PRINT 'PA' iX , 0j - 

4395 IF Ii>-0 THEN PRINT "CP -2 ,-2 ; L9 1 0 j PR -2 ,2 jLR' i I i ; " ' 

4400 IF li :0 THFN PRINT 'CP - 2 . -2 J i B I 0 j PR 0 ,2 i L B ' 1 1 ii " " 

4 405 I 1 * I i 1 

4410 NEXT I 

4415 PRINT 'PU PA 0 0' 

4420 Ii=LGT(Yninl 
4425 YI0=10 

4430 FOR 1*1 TO fy+1 
4435 Ya*Ynin* 10 ( I -1 > 

4440 Y*LPT( Ya/Ynln )«Sfy 
4445 PRINT 'PA 0 ,* i Y , " 

4450 PRINT 'CP -4 ,-.25il BIOjFR -2 „ 2 iLP " i I U ' " 

4455 Ii«ti+1 
44B0 NFXT I 
44G5 PFEP 

4470 INPUT ‘UANT U5E nFFAUI T I ABEL 5 ( 1 «Y , O-'N )3 " . I d I 
4475 IF IdlOl THEN 
4490 RFEP 

44Q5 INPUT 'E'irFR X-I.APEl.' XlahaU 
4490 BFFP 

4495 INPUT 'ENTER Y -L ARFl. “ . Y I aUe I 5 

4500 END IF 

4505 IF Ono 2 THFN 

4510 PRINT "5R l,2iPU PA 40 -U;' 

4515 PRINT 'IH(T;PR -1.6,3 IMI PR 1.2,0 PUsPR . 5 . -4 H. Bwo « PR .5,li“ 
4520 PRINT •'I.B-TiPR . 5 . - 1 i L R-^a ♦ i PR .5.M’ 

4525 IF Iun-“0 THFN 
4530 PRINT "LB I <K )' 

4535 F13;E 

4540 PRINT 'IB) < F 
4545 FMO IF 
4550 FNP IF 
4555 IF npo-2 THEN 
4560 IF hin*0 THFN 



4565 


PR INT 


"^9 1 


1 .5 ,2:PM 


PA 40 


UilBq 


( W/mSR 1,1 


.5«PR 0-5.1 II B2iSR 1 .5 ,2tPR 


0.S,- 


1 il B )■ 














4570 


El SF 














4 575 


PRINT 


"SR 1 


l.5,2iFH 


PA 34 , 


- 14 1 1 Bn 


(Btu/hrjFR 


.5. .Sil.B. iPR .5,-.5i' 



99 



45ePi PRINT “LBftiPR .5,IiSR I . I . B 1 1.82 i RR 1.5,2iPR .5,>tiLB)i* 

4588 END IF 
4590 END IF 
4595 IF Ooo=0 THEN 
4500 IF lun=0 THEN 

4R05 PRINT “SR 1.5,2iPH Pft -12,40ini 0,liLBq (U/miPR -I,0.5i 5R l.l.5ilB2iSR I 
5,2 jPR I ,.5iLB )" 

4G10 ELSE 

4G15 PRINT "SR l.5,2iPU PA -12,32501 0,liLBq (Btu/hriPR 5 . . 5 ; LB . i PR .5,.5i" 
4G20 PRINT ‘LBftiSR I.I.GsPR -1..5;LB2iPR l,.5iSR l.5,2iLB>“ 

4G25 END IF 
4G30 FND IF 
4G35 IF Opo 0 THEN 
4G40 IF lun^0 THEN 

4545 PRINT "SR l.5.2iPH PA -l2,3BiOI O.IsLBh (U/niPR -l,.5iSP I.I.5aB2iSR I. 
5.2iPR .5,.5i" 

-1650 PRINT "IB.iPR .5.0iLeK)" 

4655 ELSE 

4GG0 PRINT -SR l.5,2»PH FA -I2,29i0l 0,liLBh (Btu/hriFR > . 5 . . 5;LB . i PR .5,.5i“ 

4GG5 PRINT 'LBftiPR 'l,.5iSR l.l.5iLB25SR I.5.25PR .5,.5iL8.iPR .5,.55LBF>- 

4670 END IF 

4G75 END IF 

4680 IF Idl-0 THFN 

4685 PRINT “SR l.5.2iPU PA 50. -16 CP " i -I EN( ?( I abe I f ) /2 i "0 i LB ' i X 1 abr I f i" " 

4680 PRINT “PA -14.50 CP O , ' ; - 1 EN( fl abr 1 J )/ 2 *5 ^G i “ 0 1 0 . I iLB " i Y ) abe 1 * i " 

4G95 PRINT "CP 0.0 or 
4700 END IF 



4 705 
4710 
4715 
4720 
4725 
4730 
4735 
4740 
4745 
4?50 
4755 
47G0 
47E5 
4 770 
4775 
47B0 
47B5 
4 790 
4795 
4B00 
4805 
4B1 0 
4BI5 
4820 
4825 
4830 
4835 
4040 
4B45 
1850 
4855 
4S60 
4BG5 
4870 
4875 
4B80 
4 885 
4890 
4885 



I nn*0 



Repea t : • 

X 1 1*1 .E + G 
Xul = -I -E ^6 



ICM--0 

BEEP 

INPUT -WANT TO PLOT DATA FROM A FILE < 1 = Y . 0-N )? “ .01 
IF Ok-1 THEN 



BEEP 

INPUT “ENTER 1 HE NAME OF IHE DATA F 1 1 E ' .0_f i I r « 

ASSIGN @Fi le TO 0_fi 1 cf 

BEEP 

BEEP 

INPUT "ENTER THE BEGlNmNA RUN NUMBER", Md 
BEEP 

INPUT “ENTER THE NUMBER OF X Y PAIRS S TCREO’ ,Npa ir s 
PE£P 

INPUT “CONNECT DA 1 A UlTH LINE ( 1 *Y ,0=N )? “ . I c 1 
BEEP 

PRINTER IS I 

PRINT USING “4X, "-Select a synboi:““- 



PRINT USING “GX 1 
PRINT USING “GX "'3 
PRINT USING "BX ,“"5 
PRINT USING "GX ,"'G 
PRINT USING "GX , 
INPUT Svn 
PRINTER IS 705 

PRINT "pu or 

IF S/n-1 THEN PRINT 
IF 5vn*2 THEN PRINT 
IF Svfn = 3 rilEN PRINT 
IP f1ri ■ I THEN 
FOR n ID nifl-l ) 
ENTER i le iVa Xa 



Plii-s sign" 
Square “ 



SUr 
Circle 
RonPij a " “ " 

R Ight - s i de-up triangle"" 
Un-side-doun triangle""" 



“SM* " 
“SM» “ 
“SMo“ 



NFXT I 
END IP 

FOR I- I TO Noairs 
ENTER i«Fi I c I Ya ,Xa 



IQQ 



4900 

4305 

4910 

4915 

49:0 

49:5 

4930 

4935 

4940 

4945 

4 950 

4955 

4 9G0 

49G5 

4970 

4975 

4 980 

4 985 

4990 

4995 

5000 

5005 

5010 

5015 

50:0 

50:5 

5030 

5035 

5040 

5045 

5050 

5055 

5050 

5055 

50701 

5075 

5060 

5095 

5090 

5095 

5100 

5105 

51 10 

5115 

5i:o 

5i:5 

5130 

5135 

5140 

5145 

5150 

5155 

51G0 

51B5 

5170 

5175 

5IP0 

5165 

5190 

5195 

5:00 

5:05 

5210 

5:15 

5::0 

5::5 



IF Op 0*1 T11EH Ya«Ya/Xa 

IF Odo“2 then 

0=Ya 



Ya“Ya/Xa 

Xa»0 



END IF 

IF Xa<Xll THEN Xll-Xa 
IF Xa>XuI then Xul-Xa 
•IF Uin-I THEN 
IF Ooo<: THEN Xa«Xa*l .0 
IF OooiO THEN Ya-Ya* .1761 
IF Ooo=0 THEN Ya-Ya*.3l7 
IF O 00-2 THEN Xa-Xa*.317 
END IF 

X»LGT( Xa/Xnin )«Sf.’ 

Y-LGT<Ya/Ynin )*Sfy 
IF Svn 3 THEN PRINT '’SM" 

IF Syn- 4 THEN PRINT -SR 1. 4,2.4* 
IF IcI-0 THEN 
PRINT *PA‘ X Y . 

ELSE 



PRINT -PA- X , -PO* 
EMD IF 

IF Syn 3 THEN PRINT 
IF Syn-=4 THEN PRINT 
IF Svn=5 THEN FRINT 
IF Syn=6 THFN PRINT 
IF Svn»7 THEN PRINT 
NEXT I 
PRINI -PU- 
BEEP 

ASSIGN «Fi le TO • 
XII*XII/1 .2 



-SR 1 .: . 1 .6- 

-UC: .4 .93 ,0 ,-e .-4 ,0 ,0 ,8 4 ,01 * 
’•IJC3 0 .99, -3. -6 .-3,6 .3 .6.3.-6i" 
“UC0.5.3,99,3,-8,-6,0,3,8r’ 
-1100,-5.3,99,-3,8,6 .0,-3,-8«" 



Xul»XuI • I ,2 
GOTO 8040 
END IF 

FRINT -PU 5M- 



BEEP 

INPUT -WANT TO PLOT A POLYNOMIAL ( 1 »Y ,0=N )'? “ ,Go_on 

IF Go_cn*l THEN 

BEEP 

PRINTER IS I 

FRINT USING ’•4X ‘'*Selact line tvpe:”"' 

PRINT USING ’’6X,’“0 S^^lid Iina-“" 

PRINT USING '6X,-“1 na«5heij-"- 

PRINT L'SING "GX,“"2,, 5 Longer line - dash""" 

INPUT Ion 
PRINTER 15 705 
REEP 

INPUT “SELECT <C*LIN,1-L0G<X ,Y))- ,IIcg 
lornt = I 
CALL Pciv 

FOR X.-0 TO G.^ STEP C^/:00 



Xa"Xnin«10’Y.r 

IF Xa XII OR Xa>XuI THEN 5300 



Irn* Icn* 1 
Pu'0 

IF Inn-1 THFN Idf=Icn MOD 2 
IF lon-Z THEN Irif-Icn MOP 4 
IF Ion=3 THEN Idf=Irn MOn 0 
IF Ion-4 THEM Idf-lcn MOP IB 
IP Ion-5 THEN Idf-Icn MOD 79 
IF Idf = l 1I1FN P»j=l 
IF Ooo = 0 THEN Ya = F NPo I y < Xa ) 

IF Cno=*2 AND llog»0 THEN Va-Xa /FNPo I v ( Xa > 
IF Opo = : ANO Ilog^l rUFU ia = FNPoIy< Xa ) 
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r,2:?c 

5235 

5240 

5245 

5250 

5255 

5260 

5265 

5270 

5275 

5280 

5285 

5290 

5295 

5300 

5305 

5310 

5315 

5320 

5325 

5330 

5335 

534 0 

5345 

5350 

5355 

5360 

5365 

5370 

5375 

5380 

5395 

5390 

5395 

5400 

5405 

5410 

5415 

5420 

5425 

54 30 

5435 

54 4 0 

5445 

5450 

5455 

5460 

5465 

5470 

5475 

5480 

5485 

5490 

54 95 

5500 

5505 

5510 

5515 

5520 

5525 

5530 

5535 

5540 

5545 

5550 

5555 



IF 0oo=! THerJ Ye*FMPc] v( Xa )/Xa 

IF /a'Ymin THEM 5300 
IF Iun-1 THEM 
IF Ooo- 2 THEN Xa®<a»l.Q 
IF Ooo>0 THEN Ya-Ya» . 1761 
IF Opo=0 THEN Ya-Ya* .317 
IF Opo=2 THEN Xa»Xa« .317 
ENO IF 

Y-L6T< Ya/Ynin )*Sfv 
<»L6T( Xa.‘X'*'tn )»Sf • 

IF Y 0 THEN Y-0 

IF Y^100 THEN GOTO 5300 

IF Pu = 0 THEN PRINT “ PA " . X , Y . •PO* 

IF Pij-1 THEN PRINT “ PA" , < . Y . * PU ' 

NEXT X- 
PRINT "PU** 

ENO IF 
BEEP 

input “WANT TO QUIT ( 1 => Y , 0«N ) " . I q t 
IF Iqt=l THEN 5335 
GOTO 4715 

PRINT "FU PA 0 0 5P0" 
ciJDENO 

OEF FNHsMoct h( X ,Bco I 5d t ) 

0111 A(5) P<5 ) ,C(5 ) .0(5 ) 

DATA .20526.. 2532:.. 3 I 9048,. 55322.. 79909 I. 00253 
DATA .74515 , .72993 , . 73 1 99 , . 7 1 225 , . 68472 , . 64 I 97 
DATA .4 1092 17736 . . 25 1 42 , . 54806 . . 8 1 9! 6 , ! .0845 

PAIA . 71403, .72913,. 72565. .69RR91 665867 ,. 6 1 989 

REAO A( • ) .Of • ) ,C( • > ,D( * ) 

IF Boo 6 THEN I=8op 
IF 8oo=5 THFN 1*4 
IF Oop=10 THFN 1=5 
IF I sal"! IHEN 
H:=EXP( A( I )*-P( I T»LOG< X ) ) 

EL SE 

H3=EXPf C( I)^0( I WLOGl X ) ) 

EMP IF 
RETURN Hs 
FNENO 

DCF FnPoIv(X) 

CON /Oolv/ A( I 0 , 10 ) ,C( 10 ) .8( 4 ) Nop ,Iornl ,0po . I log 
X I»H 

PoIy*8( 0 ) 

FOR 1*1 TO Noo 

IF Ilog*l THEN X l*LCG( X ) 

PoIv=Poly + B( I )*xr I 
NEXT I 

IF IloQ“l THEN Polv=Exp(Poly ) 

RETURN Pciv 

FMEMO 

sue Poly 

OIM R( 10 ) ,5( 10 } .Sv( 12 ) ,S- ^ I 2 ) ,x>.( 100 ) ,Yv( 100 ' 

COM ''Colv/ A( 10 , 10 ) C( 10 > ,n( 4 ) N . lornt ,0oo .1 log 
CON /X-yy/ Xo( 5 ) ,Yp( 5 ) 

FOR 1*0 TO 4 
fi{ I ) = 0 
NEXT I 
8EEP 

INPUT "SELECT < 0=F I LE . 1 = KE YDOARO , 2 = PR0GRAM ) ' , In 

I n* I nx 1 

BEEP 

IMPL'T "ENTER NUM8FR Of X-Y PAIRS" .No 

IF ln«l THEN 

BEEP 

INPUT "ENTER OATA FILE NANF “ 0 file« 
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55G0 PEEP 

55B5 INPUT "LIKE TO EXCLUDE DAlA PAIRS ( I )*? " , led 

5570 IF Ied»r THEN 
5575 BEEP 

5580 INPUT "ENTER NUMBER OF PAIRS TD BE EXCLUDED" , Ipe^ 

5585 END IF 

5590 ASSIGN PFiIe TO 0_filef 
5595 ELSE 
5GO0 PEEP 

5G05 INPUT "WANT TO CREATE A DATA FILE ( l-Y ,0*N )?" ,Yes 
5GI0 IF Yes=l THEN 
5GI5 8EEP 

5G20 INPUT "6IUE A NAME FOR DATA F ILE ‘ ,D_f x I ef 
5G25 CREATE BOAT D_file* 5 
5G30 ASSIGN @FiIe TO Defile* 

5G35 END IF 
5G40 END IF 
5G45 PEEP 

5G50 INPUT "ENTER THE ORDER OF POLYNOMIAL " ,N 

5E55 FOR 1-0 TD N*2 

5GG0 Sy(I>-0 

5GG5 S>H)=0 

5G70 NEXT I 

5G75 IF led-* AND Im-I THEN 
GG80 FOR 1 = 1 TO Toe.< 

5GP5 ENTER ©File jX ,Y 
5G90 NEXT I 
5G95 END IF 
5700 FOR 1=1 TO No 
5705 IF In=l IHfN 

5710 IF Opo=: THEN ENTER ?FiIei<,Y 
5715 IF Odo'2 then ENTEP ©FiletY,X 
5770 IF I 1 00=1 then 
5775 xt=K/Y 
5730 X=LOr?(X) 

5735 Y=L06(Xt) 

5710 end if 
5745 END IF 
5750 IF Im=2 THEN 
5755 PEEP 

57G0 INPUT "ENTER NEXT X -V FAIR-^X.Y 

57G5 IF Ye5=t THEN OUTPUT ©F i 1 - i X , Y 

5770 END IF 

5775 IF In 3 THEN 

5780 X.<I)=X 

5785 Yv(I)=Y 

5790 ELSE 

5795 X*Xd<I-I) 

5B00 Y = Yp( I-l ) 

5805 END IF 
5810 R(0)=Y 

5815 Sy(0 ) = Sv(0 )fY 
5870 si 1 >--*X 
5075 S-.( I )=S^( r )<-X 

5830 FOR 1=1 TO N 
5035 R( J )-R( J-1 )*X 
5840 Sv( J ) = Sy< n+R< J ) 

5845 next J 
5850 FOP J=7 TO N»7 
5855 S( J > = S< J- 1 >‘X 
58R0 S.. ( J ) = Sn( J HS(.n 

58G5 NEXT J 
5870 NEX T I 

5875 IF Yes-1 AND ln-7 THEN 
5800 BEEP 

5095 PPINT USING ”12y,DD,“" X-Y pairs uere stored in file ‘ , 1 0A " » Np , D_f 1 1 e* 



I 
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5090 END IF 
5895 S.<(0)*No 

5900 FOR 1=0 TO N 
5905 C(I)=Sy(I) 

5910 FOR J-0 TO N 
5915 ft( I ,J )-S.k( I ►J ) 

5920 NEXT J 
59Z5 next 1 
5930 FOR 1*0 TO N-l 
5935 CrtLL Divide( I ) 

5940 CftLL Sublract< I-M ) 

5945 NFXT I 

5950 0<N)=C(N)/A(N,N> 

5955 FOR 1-0 ro M-l 
59G0 0(N- I-I >*C<N-1-I ) 

59G5 FOR J*0 TO I 

5970 D(N- I -I ) = B< N- l-I >-A(N-l -I ,N 

5975 NEXT J 

5980 0< N- I -I )*B( N- 1 -I )/A( N- 1 - I M-l -I ) 

5995 NEXT I 

5990 IPRINTER IS 701 

5995 iPPIMr S< • ) 

R0O0 (PRINTER IS 705 

6005 IF lornt*0 THEN | 

GOIO PRINT USING " 1 ZX , ‘ -EXPOflFMT CDFFF IC lENT - “ " 

6015 FOR 1=0 TO N 

6030 PRINT (15IN(3 " 1 5X .00 .5X .riD . 7DE * ; I ,B{ I ) 

6035 NEXT I 
GO30 PRINT ' - 

6035 PRINT USING “13X.‘"DATA POINT X Y Y( CALCULATED > DISCR 

EPANCY" * •• 

G04O FOR 1*1 TO Mo 
G045 Yc*B<0) 

6050 FOR J-1 TO N 
6055 Vc-Yc>0( J )*X^( I > 'J 
6060 NEXT J 
60G5 D-Yy( I ) Yc 

6070 PRINT USING • I5X .30 ,4X aiND.GDE JX )" ; I ,Xx( I ) ,Yy( I > ,Yc .0 

6075 NEXT I 

6090 END IF 

6085 ASSIGN @File TO • 

G090 SUBFND 
6095 5UD Divide<M) 

6100 CON /Colv/ A( 1 0 , 10 > ,C< 10 > ,D< 4 ) .N ,Ipr nt ,0 do .1 log 
6105 FOR I*M TO N 
G1 10 Ao-A( I ,N) 

61 15 FOR J*N TO N 

6130 A< I J )'-A( I . J ) 'Ao 

6135 NE^T J 

6130 C(n*C(I)/Ao 

6135 NEXT I 

6140 SUBFNO 

6145 SUB Sublrat t( !• > 

C150 CON /Coly/ A( 10 , 10 ) .C( 10 > „Q( 4 ) ,N , Iprnt .Ooo , I log 

6i55 FOR I*K TO N 

GIGO for J*k-i to N 

6165 A( I , J )*A(K-I J )-A{ I . J ) 

HI 70 NEXT J 

GI75 C( I )*C(K n Ct I ) 

GI90 nrxr i 

6185 SUGFND 
GI90 SUB FI in 

6195 CUN /Colv/ A( 10 . 10 ) ,C( 10 ) ,B( 4 ) ,N , lornt .Opo . 1 log 
6300 CON /X^vy/ Xx(5'.Yv{5> 

G305 PRINTER IS 70S 
6310 BEEP 
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G:15 input *GELFCT (0-h/h0% sane t ube , I *h( HF )/ h( sm ) “ , Ir t 

g::o oeep 

G::S input -UHICH T^at < 1=G.7 ,0»-2.2 >'Msal 

G230 Xt*iin*0 

G235 Xna.‘.'«10 

G240 XsteD*Z 

G245 IF lrt»0 THEN 

G250 Ynln-0 

G:-55 Yna-.«l.4 

GCGO Ysteo-.Z 

CiZGS ELSE 

G370 Ynin=0 

GZ75 Yna.<-I5 

GC60 Yst-D*5 

G295 END IF 

G290 BEEP 

G295 PRINT " INi SP I j IP 2300 ,2200 , 8300 ,GS00 « “ 

G300 PRINT "SC OJ00.0 100sTL 2 Oi" 

G305 Sfy*100/{ XnaT-Xmn ) 

G310 Sfv=l00/(Ymar'Ynin ) 

G3I5 PRINT "PU 0 ,0 PO" 

G3Z0 FOR Xa*Xnin TO Xma.^ STEP Xstep 

G325 X*< Xa-Xnin )»Sf.n 

G330 PRINT "PA'* ;X 0i XT j ' 

G335 NEXT Xa 

E340 PRINT "FA 100,0iPU«" 

G345 PRINT ‘PU PA 0,0 PD‘ 

6350 FOR Ya»Yntn TO Yna>. STEP Yst-o 
G355 Y=»< Ya-Ymin )*Sfv 

G3G0 PRINT *PA 0 “ s Y , “YT‘ 

G3G5 NEU Ya 

G370 PRINT ‘PA 0 . 1 00 TL 0 2‘ 

G375 FOR Xa = Xnin 10 Xna-. STEP Xstep 

G380 X=-< Xa-Xnin >*Sf 

G3S5 PRINT ‘FA‘ ;X . 100s X T “ 

G390 NEXT Xa 

G395 PRINT "PA 100,130 PU PA lOO.O PO“ 

G400 FOR Ya^Ynin TO Yna». STEP Y s t ep 

G405 Y*< Ya-Yn in ) »Sf y 

G410 PRINT "PD PA I00,“.Y.*YT“ 

G4iS NEXT Ya 

G420 PRINT "PA 100 , 100 FU" 

G425 PRINT 'TA 0,-2 SR I .5,2* 

6430 FOR Xa=Xnjn TO Xna\ STEP Xstep 
G435 X-< Xa-Xnin )*Sf < 

G4-10 PRINT ‘PA"!X.".0;‘ 

G445 PRINT "CP -2 I » LB " ; Xa i" * 

G450 NEXT Xa 

G455 PRINT "PU PA 0.0" 

G4G0 FOR Va=Ynin TO Ynax STEP Vsteo 
R4G5 IF ABSOa) l.E-5 THEN Ya=0 
G470 Y-( Ya-Vnln )*Sfy 

G475 PRINT "PA 0 . " ; Y , " 

G480 PRINT ‘CP -4 ,-.25iLB‘ ; Ya; “ 

G4G5 NEXT fa 

G490 X IaLieU=’‘0i I Percent" 

G495 IF Irt-0 THEN 
6500 YlabelS--"h/hOT 
6505 ELSE 

6510 V IabelS*‘h/hsnoot h* 

6515 END IP 

6520 PRINT "SR 1.5,2;PU PA 50 -10 CP " : -LEIN X labe 1 1) /2 ; " O ; I B" ; X labe 1 1 ; 
6525 PRINT ‘PA -11,50 CP 0 . “ ; - L FN( Y 1 abe 1 1 >/2 *5 ^6 ; “DI 0 , 1 iL 8 " i Y 1 abe 1 J ; 
G530 PRINT "CP 0 0" 

6535 lDr.*0 
R540 BEEP 
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6545 

6550 

6555 

6560 

6565 

6570 

6575 

6580 

G5B5 

6590 

6595 

6600 

66P5 

6610 

6615 

6620 

6625 

6630 

6635 

6640 

6645 

G650 

6655 

6660 

6665 

6670 

6675 

6680 

6685 

6693 

6695 

6700 

6705 

5710 

6715 

6720 

6725 

6720 

6735 

6740 

6745 

6750 

6755 

6'’G0 

6765 

6770 

6775 

6700 

6795 

6790 

6795 

6800 

6805 

6910 

6015 

6820 

6025 

6830 

6835 

6040 

6845 

6850 

6855 

6860 

6865 

6870 



INPUT 'WANT TO PLOT DATA FROM A FILE ( I =Y ,0-N >7 • ,01- o 



I cn«-0 

IF 01 D-1 THEN 
BEEP 

INPUT "ENTER THE NAME OF THE DATA F ILE ’ . 0_f 1 1 e* 
0EEP 

INPUT "SELECT (0-LINEAR. 1 =»LOG( X , Y > " . 1 1 og 
ASSIGN @Fi Is TO D_f i le? 

BEEP 

ENTER THE BEGINNING RUN NUMBER" 



,Md 

ENTER THE NUMBER OF X-Y PAIRS STORED * .Nqa Ir s 
ENTER DESIREO HEAT FLUX",0 



INPUT 
BEEP 
INPUT 
BEEP 
INPUT 
8EEP 

PRINTER IS 1 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
INPUT Svn 
PRINTER IS 705 
PRINT "PLI 01" 

IF Syn=I THEN PRINT 
IF Syn=2 THEN PRINT 
IF Svn=3 THEN PRINT 
Nn-4 



*4X 
"4X , 
“4X , 
"4X , 
"4X , 
‘4X , 



Select a svnbol:’"* 

1 Star 2 Plus algn*"' 

3 Circle 4 Souare*"' 

5 Rombua"’" 

6 R Ight “3 1 de-uD triangle* 

7 Ud- 5 1 de-douin triangle*’ 



”SM» " 
"snt" 
•• SMo * 



IF I Iog*l THEN Nn- 1 

IF Md.'l THFN 

FOR I't TO (Md-I ) 

ENTER PFi ieiXa .Ya 
NE/T I 
ENO IF 
OI-O 

IF I log«l THEN 0=L0G(Q ) 
FOR 1=1 TO Noairs 
ENTER ??F I lerXa ,8( • ) 
Ya=B(0) 

FOR K=1 TO Nn 



Ya= Vd+8( k )*0*K 



NEXT K 

IF I Iog=l THEN Ya>E>'P( Ya ) 
IF I Iog-0 THFN Ya»01/Yd 
IF Irt-0 THEN 
IF Xa=0 THEN 



Yo»Ya 



Ya = l 
ELSE 

Ya=Ya/To 
ENO IF 
ELSE 

Han^FNHstnoo t h( 0 ,Xa , I sa t ) 

Ya=Ya/Hsn 

ENO IF 

X - ( I-l ) = Xa 

Yv( I- I ) = Ya 

X®{ Xa-XrMn >*Sf< 

Y = ( Ya-Yrni n ) »S f v 

IF Svm 3 THEN PRINT “SM" 

IF Svm 4 THEN FPINT “SR 1 . 4 , 2 . 4 ’ 
PRINT “PA" X.y • 



IF 


Syn 3 


THEN 


PRINT 'SR 1.2,1. 


6" 


IF 


S > 4 


THEN 


PRINT "l'C2.4.99 


0 . -0 . -4 .0 .0 ,9 4 ,0 . J 


IF 


Svn = 5 


IMEN 


PRINT "UrjO 99 


.-3 .-6 .-3 .6 .3 -6 .3 .-6 



las 



6875 IF Svn>S THEN PRINT ' UC0 , 5 . 3 , 99 . 3 , - 8 . -6 0 . 3 . 8 i" 
6880 IF Syn=7 THEN PRINT “UC0 , -5 . 3 . 99 ,-3 , 8 . 6 .0 , -3 . -8 i 
6085 NEXT I 
6090 BEEP 

6895 ASSIGN gFiIe TO • 

6900 ENO IF 
6905 PRINT "PU SM* 

G9I0 BEEP 

6915 INPUT "UANT TO PLOI A POLYNONIAL (1 *Y , 0 = N )7 " , 0^ p 
6920 IF Oko'f THEN 
6925 BEEP 

6330 INPUT "SELECT ( 0-L I NEAR . 1 »L0G( X .Y ) ) ".I log 
6935 Iornt»l 
6940 CALL Poiv 

6945 FOR Xa-Xnin TO Xna.v STEP X3teo/25 
6950 Icn^IcnH 
6955 Ya=FNPoiv< >a ) 

6960 Y-( Ya-Vnin )*Sfy 

6365 X«=(Xa"Xmn >*Sf < 

6970 IF Y 0 THEN Y»0 

G975 IF Y 100 THEM GOTO 7025 

6980 Pu=0 

6985 IF Ion-1 THEN Idf=Icn MOO 2 

6990 IF Ion-2 THEN Idf=Icn MOO 4 

6995 IF Ion-3 THEN Idf-Icn MOO 8 

7000 IF Ion>4 THEN Idf-Icn MOO 16 

7005 IF Ion-5 THEN Idf-Icn MOO 32 

7010 IF Idf-I THEN Pu=I 

70)5 IF Pu=0 THEN PRINT *PA“.X Y/'PO" 

7020 IF Pu-I THEN PRINT "PA" X ,Y /'PU' 

7025 NEXT Xa 
7030 PRINT "PU" 

7035 Ipn-Ion+1 
7040 GOTO 6540 
7045 ENO IF 
7C50 BEEP 

7055 INPUT "WANT TO QUIT ( ) *Y . 0-N )7" , loul t 
7060 IF Iquit’l THEN 7070 
7065 GOTO 6540 
7070 PRINT "PU 5P0" 

7075 SUCENO 
7080 SUB Stats 
7085 PRINTER IS 70! 

7090 J-0 

7095 K=0 

7100 BEEP 

7105 INPUT "PLOT FILE TO ANALy ZE^* ,F t let 
71 10 ASSIGN @Fi le TO Filet 
7115 BEEP 

7)20 INPUT "LAST RUN No 0-OUI T >“ ,Nn 

7125 IF Nn=0 THEN 7305 

7)30 Nn-Nn-J 

7)35 Sv-0 

7140 Sy*0 

7)45 Sz-0 

7150 Sa5-0 

7!55 Sv5-0 

7160 Srs=0 

7165 FOR 1*1 TO Nn 

7170 J-Jf) 

7)75 ENTER 0Fi le;0 .T 
7180 H-O/T 
7)85 Sa-S.*.+0 

7)90 S-5-S-sfO 2 
7195 Sy-Sv+T 
7200 Sys-Sys+T 2 
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7205 Sz-52+H 

7210 Sz5-Sz3»H"2 

72IS NEXT I 
7220 Oave*Sx/Nn 
7225 Tav'e*Sy/Nn 
7230 Have-Sz/Nn 

7235 5devq-S0R( ABS( (Mn*S- 3-Sx 2 >/ ( Nn • ( Nn- 1 ) ) ) ) 

72 40 Sdevt-50R< A8S( (Nn«Sv3-Sy 2 >/(Nn«(Nn-l ))) ) 

7245 Sdevh-SOR( ABS< (Nn»Sz5-Sz 2 ) / ( Nn • ( Nn- 1 ) ) ) ) 

7250 Sh»I00*Sdevh/H3ve 
7255 So* 1 00 • Sdevo/Oave 
7250 St=l00*Sdevt/Ta^e 
7255 IF K-t THFN 7295 
7270 PRINT 

7275 PRINT USING ‘HX/'DATA F I LE : " * , I 4A ” t F i 1 
7260 PRINT 

7285 PRINT USING MIX, “•RUN Htube SdevH Odo SdevO Thetab SdevT"’ 

7290 K=l 

7295 PRINT USING "MX, DO 2( 2X .0 . 3DE . I X .30 . 2D > . 2X .DP . 3D . I X , 3D . 20" ; J .Ha v- ,Sh .Gave 
,Sq , Tave ,St 
7300 GOTO 71 15 
7305 ASSIGN ©Filel TD • 

7310 PRINTER IS 1 
7315 SUBEND 
7320 SUB Coef 

7325 con /Colv/ A< 1 0 , I 0 > C( I 0 > , B( 4 ) ,N , Iprnt ,0oo , 1 1 og 
7320 PEEP 

7335 INPUT "61UE A NAME FOR CROSS-PLDT FILE’.Cpft 
7340 CREATE BOAT Cof«,2 
7345 ASSIGN ffFl le TO Cof^ 

7350 BEEP 

7355 INPUT "SELECT ( 0-U NFAR , 1 =l 0G( X , Y > ) " . II og 
7350 PEEP 

7355 INPUT "ENTER OIL PERCENT ( - I -S TOP ) “ ,8oo 
7370 IF Bcp- 0 THEN 7390 
7375 CALL Pol, 

7360 OUTPUT @Fi leiBop ,B( • ) 

7385 GOTO 7350 

7390 ASSIGN TO • 

7395 SUBEND 
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APPEH DII C 

COBPOTBE-CONTHOLLED VAL7E PROGBAM 

The computer-controlled valve was controlled by an 
Octagon Systems SIS-2A microcontroller. The SYS-2A is a 
complete computer system on a single card and requires only 
a 5 volt supply for operation. The SYS-2A has provisions for 
4, 2.5 volt A/D data inputs and 8 high current digital 

outputs. The card was configured as shwon in Table 4 below. 



TABLE 4 

Configuration of SYS-2A I/O Channels 



Data Input Channel 
3#A10 = 8 
3#A10 = 9 



Purpose 

Valve position 

Liquid R-114 Temperature 



Data Output Channel 
S#A00 
3#A01 



Purpose 

Shut Valve Relay 
Open Valve Relay 



(0=stop, l=shut) 
(0=stop, 1=open) 



The proportional-integral-derivative control program 
below was written in NSC "tiny BASIC.” The SYS-2A microcont- 
roller operator's manual includes an appendix of tiny BASIC 
commands. Remarks are provided below that were not in the 
actual program (to minimize execution time) for clarity. The 
program was stored in ESPROM added to the card after program 
debugging was complete. 

10 Let A=15:Let B=5: Let C=lO;Let F=0:Let 1=0 Initialize 
15 Let 3#A00=0:Let S#A01=0 Stop Valve 

20 Print”Set Constants (1=Y, 0 = N) ?”:Input Q 
30 If(Q = 0)then GOTO 60 
40 Print”A,B,C,G?”:Input H,J,K,L 
50 Let A=H:Let B=J:Let C=K: Let G=L 

60 Print”Input Desired Temp (-1=End) ”:Input R Enter Temp 
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70 Let S#A10=9 ;Delay 2:Let M=1 6*o)#A1 1+S#A12/16 Read Channel 9 
80 Let 3#A 10=8: Delay 2:Let P = 1 6*3 #A 11 + 3# A 1 2/16 Read Channel 8 
90 Print"Temp=*’,M,"Pos=",P 

100 If{R<0)then GOTO 800 Exit Program 

200 Let 3#A10=8;Delay 2;Let P=1 6*3#A 1 1 + 3#A 1 2/1 6 Read Channel ’s 
210 Let 3=0 Ave 10 Readings 

220 For N=1 to 10 Step 1 



230 
240 
250 
260 
262 
264 
266 
26 8 
270 
300 
310 
400 
410 
420 
430 
440 
450 
460 
465 
470 
480 
490 
500 
510 

52 0 

53 0 
540 
60 0 



Let 3#AlO=9:Delay 2:Let H=16*3#A 1 1+3 #A 12/1 6 Read Temp 



Let E=M-R:Let S=S+E 
Next N 

let E=S:Let D=S-F:Let F=E 
If (E<0) then 1=1- 1 
If(E>0)then 1=1+1 
If (K-31000) then I=-31000 
If (I>31000) then 1=31000 
Print"M='',H,»E= '’D=», D, ”I=», I 
If (E<- 1000) then GOTO 600 
If (E>2000) then GOTO 700 
Let V= (E/A)+ (B*D) + (I/C) 

Print"V=”, V, ’»ET=”,E/A,"DT=»,B*D, "IT=»,I/C 
If(V>0)then GOTO 500 
Let V=V*{-1)*G 



Error; average 

Error ; Derivative 
Integral 
Integral 
Binary Limit 
Binary Limit 



Control Band 
Control Band 
PID Value 
Debug Tool 
Above Desired 
Shut Faster 
Check 

Delay Function limited 1-1040 entry 
Slow Valve near end of travel 



If(V<1)then V=1 
If (V>1040) then V=0 
If (P<150) then V=15 
If(P<50)then V=10 

If(P<20)then GOTO 200 End of Travel Limit 

Let 3#A00=1:Delay (V) :Let 3#A00=0 Shut Variable Time 

GOTO 200 

If (P<2500) then GOTO 200 End of Travel Limit 
If(V<1)then V=1 Check 

If (V>1040) then V=0 Delay Function limited 1-1040 entry 
Let 3#A01 = 1 ; Delay (V) ; Let 3#A01 = 0 Open Variable Time 

GOTO 200 

If(P<20)then GOTO 200 End of Travel Limit 
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Out of Bound Warning 
Max valve shutting speed 
Slow near end of travel limit 



610 Print”600 loop" 

62 0 Let V=0 

630 If (P<250) then V=50 
640 If (P<150) then V=15 
645 If(P<50)then V=10 
650 let a#A00= 1 ; Delay {V) : Let 3#A00=0 Shut Valve Fully 

660 let a#AlO=8:Delay 2:Let P=16*a#A11+3#Al2/16 Read Channel 8 
670 GOTO 600 

700 If (P>2500) then GOTO 200 End of Travel limit 

710 Print"700 loop" Out of Bound Warning 

720 let a#A0 1= 1 ; Delay 0;Let S#A01=1 Open Valve Fully 

730 let a#AlO=8:Delay 2:Let P=16*3#Al1+a#Al2/16 Read Channel 8 

740 GOTO 700 

800 END 
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APPENDIX D 



EXiaPLES OP REPBESEHTATI7E DATA RONS 

The two printouts below are samples of the output of 
DRP2. The first printout (ffH79) is for the smooth tube at 
high heat flux with 0 percent oil. The second printout 
(HF101) is for the High Flux tube under similiar conditions. 



Month, date and time :03: 25 : 1 0 : 1 8 

NOTE: Program name : 0RP2 
Disk number » 02 
Old file name: UH75 

This data set ta^ en on : 02:10:12:50:15 
Tube Tyoe: 4 



Data Set 


Number 


- 1 Bull Oil r, = 


0.0 


02; 


10:13:16:18 


TC No: 


1 


2 3 4 


5 


6 


7 9 


T emo : 


IG.42 


14. G9 14.01 15.00 


14. G7 


14.91 


16.79 14.74 


Twa 


Tl i od 


Tliod2 Tvaor Psat 


T sumo 






15.23 


-2.26 


-2.20 -.63 -6.31 


-17.7 






Thetab 


Htube 


Odo 








1G.935 


5.704E^03 9.GG0E+04 








Oata Set 


Number 


- 2 8ulK Oil X - 


0.0 


02: 


10: 13: 19:35 


TC No: 


1 


2 3 4 


5 


6 


7 8 


T emo : 


16.40 


14.60 14.54 14.93 


14.52 


14.74 


16.69 14.67 


Tua 


T 1 1 Qd 


T I iod2 Tvaor Psat 


T sumo 






15.15 


-2.27 


-2.10 -1.44 -G.32 


-17.1 






Thetab 


Htube 


Odo 








16.866 


5.745e-»-03 g.689£-^04 








Oata Set 


Number 


- 3 Bulk Oil X - 


0.0 


02: 


10:13:19:56 


TC No: 


1 


2 3 4 


5 


6 


7 9 


T emo : 


16.48 


14.56 14.52 14.95 


14.50 


14.67 


16.62 14.67 


Tua 


T 1 1 od 


Tliod2 Tvaor Psat 


T sumo 






15.12 


-2.27 


-2.13 -1.51 -6.32 


-17.0 






Thetab 


Htuoe 


Odo 








16.840 


5.7G0E4-03 9.701E>04 








NOTE: 03 


X-Y oairs were stored in plot data 


file 


PUH79 



DATA FILE:PUH79 

RUN Htube SdevH Odo SdevO Thetab SdevT 

3 5.73GE+03 .51 9.G83E404 .22 1G.880 .29 
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Month, date and tine : 03 : 25 : ) 0 : 03 : 2G 

NOTE: Progran nane : 0RP2 
Oisk nuntaer = 02 
Old file nane: HF I 0 1 

This data set taken on : 02 : I S ; 1 9 : 06 : 52 
Tube Type: 5 



Data Set 


Nunber 


* 1 Bulk Oil X * 


0.0 


02:16: 19:31 :02 


TC No: 


1 


2 3 


4 


5 


6 7 8 


T eno : 


5.52 


2.63 3,76 


4.57 


3.93 


4.52 7.31 4.17 


Tua 


T 1 iqd 


Tliad2‘ Tvapr 


Psdt 


Tsuna 




4.55 


-2.14 


-1.87 .37 


-6.20 


-15.7 




Thetab 


Ht ube 


Odo 








3.370 


2.801E+04 9.440E+04 








Data Set 


Nunber 


“ 2 Bulk Oil X *= 


0.0 


02: 16: 19:31:16 


TC No: 


1 


2 3 


4 


5 


6 7 8 


Tena : 


5.54 


2,65 3.76 


4.57 


3.93 


4,53 7,31 4.18 


Twa 


TUqd 


T 1 iqd2 Tvapr 


Psat 


T sunp 




4.56 


-2.14 


-1.90 .32 


-6.20 


-15.6 




Thetab 


Htube 


Qda 








3.379 


2.792E+04 9.433E+04 








Data Set 


Nunber 


- 3 Sulk, 01 1 X - 


0,0 


02:16: 19:31 :27 


TC No: 


1 


2 3 


4 


5 


6 7 9 


Tena : 


5.55 


2.63 3.76 


4.57 


3.93 


4.53 7.31 4.17 


T Wd 


Tliad 


T 1 iad2 T vaor 


Psat 


T suno 




4.56 


-2.14 


-1,91 ,26 


-6.20 


-15.5 




Thetab 


Htube 


Oda 









3.369 2.807E+04 9.456E+04 

NOTE: 03 X-Y pairs were stored in plot data file PHFI0I 



DATA FILE:PHFI01 

RUN Hlube SdevH Qda SdevQ Thetab SdevT 

3 2.800E+04 .27 9.443E+04 .12 3.373 .16 
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APPENDIX E 
ONCERTAINTY ANALYSIS 

The uncertainty of the heat-transfer coefficient at 37 
kW/m 2 and 5 kW/m^ of runs HH79 and HF101 are analyzed below. 
The analysis is based on the Kline-McClintock [Be£. 24] 
method of uncertainty analysis. 

The heat-transfer coefficient is: 



h 



qc 

^wo ” ^sat 



and 



(E.1) 



I 



T - T = fwl - (02/0^) _ ^ 

wo sat 2 TT k L sat 



(E.2) 



where 

h 

_^c 
*^wo 
^ sat 
^wi 
Qc 
D 2 
Dl 
k 
L 

let 

F 



heat- transfer coefficient 
heat flux corrected for end losses 
average outer wall temperature 
saturation temperature 
average inner wall temperature 
heat input corrected for end losses 
outer wall diameter of tube 
inner wall diameter of tube 
thermal conductivity of wall 
length of heated surface 



Qc (1)2 /Dx) 



2 ir k L 



(E . 3 } 



According to Kline and HcClintock, the uncertainty in the 
heat-transf er coefficient is: 




fcy neglecting the error from the logarithmic term, because 
it is small compared to the other terms, the uncertainty of 
the Fourier term (F) can be estimated as: 




Oc = qc TT r>2 L 

which has an uncertainty of 



(E.5) 



(E.6) 




(E.7) 



Table 5 lists the various terms of equations (E.1) through 
(E.7) assuming all 8 wall thermocouples are used to get the 
average wall temperature f • . The data-reduct ion program 
(DRP2) used this method to calculate the heat- transfer 
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TABLE 5 

Oncertainty Analysis Terms Osing 8 Thermocouples 



File 

Heat Flux 

60 c 



6k 

~k 

6L 

“l 



WH79 

37 kW/ni2 



0.007 



HH79 HF101 

5 kW/m2 37 kV/m^ 



0.15 



0.0005 



F term (OC) 0.201 



0-009 

0.15 

0.0005 

0.024 



0.006 

0.33 

0-0005 

1.316 



HF1 01 
5 JcH/m2 



0. 007 
0.33 
0. 0005 
0. 170 



S? (”C) 


0.030 


0.004 


0.434 


0.056 




11.46 


5.63 


0.79 


-1.43 




-2.21 


-2.23 


-2.21 


-2. 18 


«T„1 


0.032 


0-073 


0.290 


0. 126 


^wo ^sat 


6F 


0.002 


0.0005 


0. 259 


0. 097 


T - T 

^wo ^sat 


^sat 


0-003 


0.008 


0.015 


0.019 


T - T 

•■wo ^sat 




0.003 


0.007 


0.008 


0- 003 


6h 


h 


0.032 


0.074 


0.389 


0. 160 


h (W/m2 K) 


2660 


540 


22300 


8400 
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coefficient. Table 6 lists the same terms assuming only the 
center 4 wall thermocouples (2, 3, 5, and 6) are used to 
calculate the heat-transfer coefficient. Comparing Tables 5 
and 6 shows the effect of the axial and circumferential wall 
temperature distributions on the uncertainty of the 
heat-transfer coefficient. 

The constraining error of the smooth tube is the uncer- 
tainty in the wall temperature. Removing the effect of the 
axial wall temperatue distribution reduces the uncertainty 
by 1 percent. The large wall superheat of the smooth tube 
contributes to the small magnitude of the uncertainty terms. 
The constraining error of the High Flux tube is also the 
uncertainty in the wall temperature, but the uncertainty of 
the thermal conductivity (part of the F term) is about the 
same magnitude an results in the larger overall uncertainty 
of the High Flux data. The small wall superheats of the 
High Flux tube also amplify the magnitudes of uncertainty 
terms. Removing the effect of the axial wall temperature 
distribution makes the uncertainty in the thermal conduc- 
tivity of the copper- nickel High Flux tube the constraining 
uncertainty. The axial temperature distribution is respon- 
sible for about 25 percent of the wall uncertainty term, but 
again the combined effect of the uncertainty of the thermal 
conductivity, and low wall superheats, does not make it 
wholely responsible for the large uncertainty of the High 
Flux data. More accurate data could be obtained on the High 
Flux surface by using a solid copper tube without a large 
uncertainty in the wall resistance or an axial wall 
temperature variation. 
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TABLE 6 



Uncertainty 


Analysis 


Terms Using 


Center 4 


Thermocouples 


File 

Heat Flux 


WH79 

37 kW/m2 


WH79 
5 kW/m2 


HF101 
37 kH/m 


HF101 

2 5 kH/m2 


50c 

Qc 


0.007 


0.009 


0.006 


0. 007 


6k 

T 


0.15 


0. 15 


0.33 


0.33 


6L 


0.0005 


0.0005 


0.0005 


0.0005 


F term (°C) 


0.201 


0.024 


1.316 


0. 170 


6F (og 


0.030 


0.0 04 


0.434 


0. 056 


f (OC) 

wi 


11.34 


5.80 


0. 44 


-1.49 


T (OC) 

sat 


-2.21 


-2.23 


-2.21 


-2. 18 


5Twi 
"" "^sat 


0.015 


0.063 


0- 21 3 


0. 097 


6F 

“ "^sat 


0.002 


0. 0004 


0. 325 


0. 108 


^sat 

"^wo " ^sat 


0.003 


0.008 


0.019 


0. 021 


6qc 

qc 


0.003 


0.007 


0.008 


0.003 


6h 

“h 


0.0 16 


0.064 


0.389 


0. 146 


h (W/m2 K) 


2690 


530 


27800 


9280 
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